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FOREWORD 


This  final  report  describes  the  research  performed  during  the  past  year  on 
adaptive  antenna  techniques.  The  research  included  new  work  on  the  topic  of 
signal  cancellation  in  adaptive  arrays  as  well  as  a  further  investigation  of 
power-separation  techniques  that  wras  begun  under  a  prior  contract  (N00019-- 
78-C-0276)  with  the  Naval  Air  Systems  Command.  The  report  has  been  divided 
into  two  parts  to  reflect  the  different  thrusts  of  the  new  and  the  continuing 
investigations. 

Part  I  of  the  report  is  concerned  with  the  phenomenon  of  signal  cancella¬ 
tion  that  occurs  in  standard  adaptive  beamformers  when  the  environment 
demands  rapid  adaptation.  The  signal-cancellation  effect  that  results  from 
jammer/signal  interaction  with  the  adaptive  beamformer  is  described  and  is 
likened  to  the  notch-filtering  effect  that  occurs  in  much  simple  adaptive  noise 
cancelling  systems.  Two  different  techniques  that  are  currently  being  pursued 
to  alleviate  signal  cancellation  are  described,  and  simulation  results  are 
presented  that  indicate  significant  peformance  improvement  over  conventional 
beamformers. 


Part  II  is  based  on  a  doctoral  dissertation  by  T.  Saxe  and  represents  a 
refinement  of  earlier  work  on  adaptive  power  separators.  This  part  discusses 
the  design  of  adaptive  power  separators  for  discriminating  between  powerful 
jammer  signals  and  a  weaker  desired  signal.  Three  different  single-channel, 
lo'.v-power-pass  adaptive  separators  are  analyzed.  Each  structure  consists  of  a 
power-spcctrum  estimator,  a  filter  designer  that  analyzes  the  power  spectrum 
and  designs  filters  that  reject  frequencies  wrhere  the  spectral  estimate  exceeds 
a  selected,  threshold,  and  a  filter  that  actually  processes  the  input  to  reject 


powerful  spectral  components. 
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I.  INTRODUCTION 


Adaptive  antennas  have  been  under  development  in  various  forms  during 
the  past  two  decades  or  so.  Although  adaptive  antennas  have  only  been  used  in 
small  numbers  thus  far,  they  have  proven  themselves  capable  of  rejecting  jam¬ 
ming  signals  to  an  extent  that  is  unprecedented.  Most  high  performance  radar 
and  communications  systems  being  designed  to  work  in  jamming  environments 
currently  incorporate  adaptive  antennas.  The  development  of  spread  spectrum 
techniques  simultaneously  with  adaptive  antennas  provides  a  formidable  set  of 
technologies  for  jam  resistant  systems.  These  technologies  are  compatible  and 
are  frequently  used  in  the  same  system.  An  adaptive  antenna  is  relied  upon  to 
attenuate  strong  jamming  signals  as  they  appear  at  the  receiver  ’’front  end.” 
Spread  spectrum  techniques  are  used  to  neutralize  large  numbers  of  weak  jam¬ 
mers  that  may  not  be  eliminated  totally  by  the  adaptive  antenna. 

A 

The  question  arises,  are  adaptive  antennas  susceptible  to  certain  jam¬ 
mers?  Or,  stated  differently,  can  jammers  be  devised  specifically  for  use 
against  adaptive  antennas?  Under  certain  circumstances,  the  answer  to  this 
question  is  yes.  This  paper  is  concerned  with  jamming  signals  that  could  defeat 
or  partially  defeat  known  adaptive  antenna  algorithms.  The  existence  of  jam¬ 
mers  that  could  trouble  known  adaptive  arrays  motivates  the  development  of 
new  adaptive  signal  processing  and  array  processing  algorithms,  two  of  which 
are  proposed  herein. 

The  goals  of  this  paper  are  threefold: 

(a)  Examine  the  signal  cancellation  phenomenon  in  adaptive  beamform- 
ers. 

(b)  Formulate  approaches  toward  the  elimination  of  the  adaptive  signal 
cancellation  phenomenon,  based  on  the  work  of  K.  Duvall. 
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(c)  Introduce  spatial  dither  algorithms  for  the  purposes  of  combating  sig¬ 
nal  cancellation  and  modulating  "smart"  jammer  signals  at  the  receiv¬ 
ing  array. 


II.  SIGNAL  CANCELLATION  JAMMING 

Any  adaptive  beamformer,  either  the  Howells-Applebaum,  sidelobe  canceller 
[1.2],  Widrow’s  pilot  signal  beamformer  [3],  Griffith’s  beamformer  [4],  Frost’s 
beamformer  [5],  Zahm’s  beamformer  [6]  or  combinations  and  variations  of 
these,  is  susceptible  to  attack  by  a  simple  jammer  which  may  be  band-pass 
noise,  or  a  sinusoid,  or  a  sum  of  sinusoids  suitably  spaced  in  frequency.  The 
interaction  of  such  jammers  with  the  desired  signal  in  these  adaptive  algo¬ 
rithms  can  cause  cancellation  of  signal  components,  even  when  these  adaptive 
beamformers  are  working  perfectly. 

To  understand  how  this  comes  about,  consider  the  Frost  beamformer  which 
functions  in  the  following  manner.  A  beam  is  formed  toward  a  user  selected 
"look  direction."  The  receiving  sensitivity  in  this  direction  is  then  constrained. 
A  typical  constraint  is  one  that  forces  the  array  to  have  a  unit  gain  magnitude 
and  zero  phase  over  a  selected  passband  of  frequencies  in  the  look  direction. 
The  beamformer  is  adapted  (its  weights  are  varied)  to  minimize  its  output 
power,  subject  to  the  constraint  which  sustains  the  beam  in  the  look  direction. 
Adaptation  subject  to  the  constraint  causes  the  array  to  accept  a  signal  with 
gain  one  if  this  signal  arrives  from  the  look  direction,  and  causes  any  other  sig¬ 
nals,  jammer  signals  for  example,  to  be  rejected  as  well  as  possible  (in  the 
minimum  total  power  sense)  so  long  as  they  do  not  arrive  from  the  look  direc¬ 
tion.  Other  adaptive  beamformers  behave  more  or  less  like  the  Frost  beam- 
former  except  that  the  Frost  algorithm  imposes  a  "hard"  constraint  on  the  sig¬ 
nal  gain  in  the  look  direction.  The  Widrow  and  Griffiths  beamformers  create 


"soft"  constraints  in  this  direction.  The  Howells-Applebaum  and  Zahm 


bcainformers  apply  soft  constraints  omnidirectionally  rather  than  along  a  look 
direction. 

Suppose  that  the  Frost  beamformer  has  a  sinusoidal  input  signal  arriving 
from  the  look  direction.  This  signal  should  appear  at  the  beamformer  output, 
going  through  a  unit  gain.  Now  suppose  a  jammer  is  turned  on,  a  very  strong 
sinusoidal  jammer  at  the  same  frequency  as  the  signal  and  arriving  ofT  the  look 
direction.  The  jamming  sinusoid  would  normally  be  rejected  by  the  adaptive 
beamformer,  if  the  signal  were  not  present.  But  with  signal  present,  minimiz- 
ing  the  total  output  power  will  cause  the  jammer  to  be  admitted  with  just  the 
right  magnitude  and  phase  to  cancel  the  sinusoidat  signal.  Thus,  the  signal 
sinusoid  is  admitted  with  a  gain  of  one.  On  the  other  hand,  just  a  trickle  of  the 
powerful  jammer  sinusoid  is  admitted  to  perfectly  cancel  the  signal  sinusoid 
and  produce  a  net  output  of  zero.  The  output  power  is  minimized  and  the  con¬ 
straint  is  preserved,  as  it  should  be  with  a  perfectly  working  Frost  beamformer. 
But  the  result  is  loss  of  the  signal.  This  amounts  to  jamming  by  signal  cancel¬ 
lation.  rather  than  jamming  by  overwhelming  the  signal  with  interference. 

If  the  input  signal  in  the  look  direction  is  broadband  rather  than  sinusoidal 

and  the  jammer  is  sinusoidal,  the  adaptive  algorithm  will  modulate  the 

sinusoidal  jammer  in  such  a  manner  that  it  will  cancel  some  signal  components 

at  the  jammer  frequency  and  at  neighboring  frequencies.  If  the  jammer  signal 

contains  a  sum  of  sinusoids  at  spaced  frequencies  within  the  passbands,  the 

output  signal  spectrum  will  be  notched  at  each  of  the  jammer  frequencies. 

This  phenomenon  could  be  troublesome  for  bandpass  and  spread  spectrum  com- 
s 

munications. 

This  signal  cancellation  phenomenon  has  been  observed  and  analyzed  in 
the  context  of  simple  adaptive  noise  cancelling  systems,  much  simpler  systems 
than  adaptive  beamformers.  A  brief  discussion  and  analysis  of  adaptive  noise 


cancelling  follows. 


HI.  ADAPTIVE  NOISE  CANCELLING 

An  adaptive  noise  canceller  is  shown  in  Fig.  1.  In  the  terminology  of  that 
field,  the  "primary  input"  contains  a  useful  signal  s,  plus  interference  Tie.  The 
"reference  input”  is  separately  obtained  in  practical  systems.  It  contains 
interference  nj,  related  to  that  of  the  primary  input.  The  relationship  between 
the  two  interferences  is  generally  unknown  a  priori.  The  adaptive  filter  has  the 
job  of  shaping  the  reference  interference  to  replicate  (in  the  least  squared 
error  sense)  the  primary  interference  so  that  subtraction  will  remove  the 
interference  from  the  primary  input  and  thereby  deliver  a  much  more  useful 
output.  It  has  been  shown  in  [7]  that  an  adaptive  filter  minimizing  mean 
squared  error,  minimizes  output  power  of  the  system  of  Fig.  1.  This  causes  the 
system  output  to  be  a  best  least  squares  estimate  of  the  useful  signal  s.  The 
Howells-Applebaum  sidelobe  canceller  is  more  complicated  than  this  in  several 
ways  (useful  signals  and  jammer  signals  appear  at  both  primary  and  reference 
inputs;  spatial,  i.e.  array,  processing  is  also  involved),  but  works  basically  on 
the  cancelling  principle  described  above. 

If  the  reference  input  is  a  sinusoid,  as  shown  in  Fig.  2,  then  the  signal  flow- 
path  from  primary  input  to  the  output  behaves  like  a  sharp,  linear,  time- 
invariant  notch  filter.  When  first  discovered,  this  was  a  surprise  because  the 
adaptive  filter  itself  is  intrinsically  nonlinear  and  time  variable.  An  analysis,  by 
John  Glover  of  the  notch  filter  effect  was  presented  in  a  1975  Proceedings  of  the 
1EIJE  paper  by  Widrow,  et  al,  on  the  subject  of  adaptive  noise  cancelling  [7].  A 
more  detailed  analysis  is  contained  in  Glover’s  Ph.D.  thesis  entitled  "Adaptive 
Noise  Cancelling  of  Sinusoidal  Interferences",  Department  of  Electrical 
Engineering,  Stanford  University,  1975  [8].  A  published  work  based  on  Glover’s 
thesis  has  since  appeared  [9].  This  work  treats  both  single  and  multiple  notch 
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cases.  Analysis  of  the  simplest  case,  a  single  notch  created  by  a  two  weight 
adaptive  filter,  is  presented  next. 

IV.  AN  ADAPTIVE  NOTCH  FILTER 

In  this  section,  an  analysis  of  the  notch  filter  effect  of  the  adaptive  noise 
canceller,  is  presented.  This  analysis  deals  with  the  formation  of  a  notch  at  a 
single  frequency.  Analytical  and  experimental  results  show,  however,  that  if 
more  than  one  frequency  is  present  in  the  reference  input,  a  notch  for  each 
frequency  will  be  formed. 

Figure  2  shows  an  adaptive  noise  canceller  with  two  adaptive  weights.  The 
primary  input  is  assumed  to  be  an  arbitrary  signal-stochastic,  deterministic, 
periodic,  transient,  etc.  The  reference  input  is  assumed  to  be  a  pure  cosine 
wave  C  cos(cj„< +<p).  The  primary  and  reference  inputs  are  sampled  at  the  fre¬ 
quency  0  =  2 i\/T  rad/s.  The  reference  input  is  sampled  directly,  giving  xj^,  and 
after  undergoing  a  90°  phase  shift,  it  is  sampled,  giving  Xzj-  The  samplers  are 
synchronous  and  strobe  at  t  =0,  T,  2T,  etc. 

A  transfer  function  for  the  noise  canceller  of  Fig.  2  may  be  obtained  by 
analyzing  signal  propagation  from  the  primary  input  to  the  system  output.* 
Weight  updating  in  the  system  is  carried  out  according  to  the  LMS  algorithm 
[10.11]: 

Wlj+i  -  wif  + 

iuz*+i  =  wS}  +  Ztir.jXy  ,  (1) 

* 

where  the  subscripts  indicate  the  time  index  and  /jl  is  a  constant  controlling  the 
rate  of  adaptation.  Refering  to  Fig.  3,  the  sampled  reference  inputs  are 

*.'t  is  not  obvious,  from  inspection  of  Fig.  2,  that  a  transfer  function  for  this  propagation 
p.-.th  in  fact  exists  Its  existence  Is  shown,  however,  by  the  subsequent  analysis. 


Xij  =  C  cos(cj0jT +d») . 


(2) 


and 

1 2j  =  c  sin(o0jT +p) .  (3) 

The  first  step  in  the  analysis  is  to  obtain  the  isolated  impulse  response 
from  the  error  tj,  point  C,  to  the  filter  output,  point  G.  with  the  feedback  loop 
from  point  G  to  point  B  broken.  Let  an  impulse  of  amplitude  a  be  applied  at 
point  C  at  discrete  time  j  =  k\  that  is, 

tj  =  a6(j-k)  (4) 

The  6{j-k)  is  a  Kronecker  delta  function,  defined  as 


<5(j-Jfc)  =  for  j=k 

J  0,  otherwise  . 

The  impulse  causes  a  response  at  point  D  of 


(5) 


_  aC  cos(v0kT+if>),  for  j  =  k 
-  * 

0,  otherwise , 


(6) 


which  is  the  input  impulse  scaled  in  amplitude  by  the  instantaneous  value  of 
xlj  at  j  —  fc*  The  signal  flow  path  from  point  D  to  point  E  is  that  of  a  digital 
integrator  with  transfer  function  2/V(z~l)  and  impulse  response  2fj.u(j- 1). 
where  u(j)  is  the  discrete  unit  step  function 


u 


O')  = 


0. 

1, 


for  j  <  0 
for  j  £  0. 


(7) 


Convolving  2fiu(J-l)  with  tjiXj  yields  a  response  at  point  E  of 

'■  w\j  =  2fj.aC  cos(u0kT+<p)  ,  (B) 

where  j  £  k  +  1.  When  the  scaled  and  delayed  step  function  is  multiplied  by 
x\j,  the  response  at  point  F  is  obtained  as 


•yij  =  2^aC*  cos(cj0jT+^)  cos(w0fc7’+p)  . 


(9) 


where  j  J  fc  +  1.  The  corresponding  response  at  point  J,  obtained  in  a  similar 
manner,  is 


y zj  -  2fiaC*  sin(u0jT  +  p)  sin(uekT+y)  ,  (10) 

where  j  k  +  1.  Combining  (9)  and  (10)  yields  the  response  at  the  filter  out¬ 
put,  point  G: 


Vzj  ~  2/iaC*  cos(a0  T(j-k  )) 

~  2fiaC*u(j-k-l)  cos(u0T(j-k))  .  (11) 

* 

Note  that  (11)  is  a  function  only  of  O'— fc)  and  is  thus  a  time  invariant  impulse 
response,  proportional  to  the  input  impulse. 

A  linear  transfer  function  for  the  noise  canceller  may  now  be  derived  in 
the  following  manner.  If  the  time  k  is  set  equal  to  zero,  the  unit  impulse 
response  of  the  linear  time-invariant  signal-flow  path  from  point  C  to  point  G  is 

Vj  =  2jiCzu(j-l)  cos(q0jT)  ,  (12) 


and  the  transfer  function  of  this  path  is 

z(z  -  cos(y07’) 


C(z)  =  2nCz 


zz  -  2zcosu07’  +  1 
2fiCe(z  cos  w0T-\) 
z2  -  Zz  cos  u0 T  +  1 


-  1 


(13) 


This  function  can  be  expressed  in  terms  of  a  radian  sampling  frequency 
f]  =  2tt/T  as 


2fiCz[z  cos(2TrcJon_l)  -  1  ) 
z2  -  2z  cos(2ncjon-1)  +  1 


(14) 


If  the  feedback  loop  from  point  G  to  point  B  is  now  closed,  the  transfer  function 
K(z)  from  the  primary  input,  point  A.  to  the  noise  canceller  output,  point  C,  can 
be  obtained  from  the  feedback  formula: 


M(z) 


z2  -  2z  cos(2n«,n“1)  +  1 
z2  -  2(1  -fxCz)z  cos(2n&3en_I)  +  1  -  2fiCz 


(15) 
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Equation  (15)  shows  that  the  noise  canceller  with  a  cosine  reference  input 
has  the  properties  of  a  notch  filter  at  the  reference  frequency  o0  along  the  sig¬ 
nal  flow  path  from  primary  input  to  output.  The  zeros  of  the  transfer  function 
are  located  in  the  Z  plane  at 

z  =  exp  (±a2rroan_1)  (16) 

and  are  precisely  on  the  unit  circle  at  angles  of  ±2nuBQ~t  rad.  The  poles  are 
located  at 

z  =  (1-/zC2)cos(27tw0-1)  ±i[(l-2 fj.Cz)  -  (l-/zC2)cos2(2ncJon-,)],'/2  .  (17) 


The  poles  are  inside  the  unit  circle  at  a  radial  distance  (l  —  S/zC2)!/2  ,  approxi¬ 
mately  equal  to  1  -  iiCz  ,  from  the  origin  at  angles  of 

±arc  cos[(l-/zC2)(l-2p;C2)~1/z  cos^rr/zO'1)]  . 

For  slow  adaptation  (that  is,  small  values  of  //C2)  these  angles  depend  on  the 
factor 


1-//C2 
(1—2  mC2)1/2 


1  -  2 piC*  4-  /x2C4 
(  1  -  2 fiCz 

(i  -  fj.zc*  +  •  •  •  ys* 


(18) 


which  differs  only  slightly  from  a  value  of  one.  The  result  is  that,  in  practical 
instances,  the  angles  of  the  poles  are  almost  identical  to  those  of  the  zeros. 

The  location  of  the  poles  and  zeros  and  the  magnitude  of  the  transfer  func¬ 
tion  in  terms  of  frequency  are  shown  in  Fig.  4.  Since  the  zeros  lie  on  the  unit 
circle,  the  depth  of  the  notch  in  the  transfer  function  is  infinite  at  the  fre¬ 
quency  u  =  wc.  The  sharpness  of  the  notch  is  detemined  by  the  closeness  of 
the  poles  to  the  zeros.  Corresponding  poles  and  zeros  are  separated  by  a  dis- 
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tance  approximately  equal  to  /zC2.  The  arc  length  along  the  unit  circle  (cen¬ 
tered  at  the  position  of  a  zero)  spanning  the  distance  between  half-power  points 
is  approximately  2/zC8.  This  length  corresponds  to  a  notch  bandwidth  of 

(BW)  =  fxCzQ/n 

=  ZfxC*F  Hz  .  (19) 

where  F  is  the  sampling  frequency  in  Hz.  The  Q  of  the  notch  is  determined  by 
the  ratio  of  the  center  frequency  to  the  bandwidth. 


fiC20' 


The  time  constant  of  the  mean  square  error  "learning  curve"  for  the  LMS 


algorithm  has  been  shown  to  be  [6,9] 


4 ju  trace  R 


iterations. 


where  R  is  the  covariance  matrix  of  the  inputs  to  the  weights,  and  n  is  the 
number  of  weights.  Formula  (21)  applies  when  the  eigenvalues  are  all  equal. 
This  is  the  case  for  the  system  of  Fig.  2.  Multiplying  by  the  sampling  period  T, 
the  time  constant  is  expressed  in  seconds  of  real  time  as 

T—  =  5oc-  (22> 

For  the  two-weight  adaptive  filter  of  Fig.  2, 

trace  R  =  +  -i<72 

2  2 

=  C2 .  (23) 

This  is  the  sum  of  the  power  into  the  weights.  Combining  equations  (23).  (22). 


and  (19)  yields 


(BW)  =  Hz. 

TmH 


Thus,  the  bandwidth  of  the  notch  is  the  reciprocal  of  the  time  constant  of  the 
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t.'ince  approximately  equal  to  /zC2.  The  arc  length  along  the  unit  circle  (cen¬ 
tered  at  the  position  of  a  zero)  spanning  the  distance  between  half-power  points 
is  approximately  2/zC®.  This  length  corresponds  to  a  notch  bandwidth  of 

(B  W)  =  ftCzQ/-n 

=  2{iC*F  Hz  .  (19) 

where  F  is  the  sampling  frequency  in  Hz.  The  Q  of  the  notch  is  determined  by 
the  ratio  of  the  center  frequency  to  the  bandwidth. 


CJ07I 

~jlczn' 


(20) 


The  time  constant  of  the  mean  square  error  "learning  curve"  for  the  LMS 
algorithm  has  been  shown  to  be  [6.9] 


Tmse  ~  — T~ - iterations,  (21) 

mse  4/z  trace  R 

where  R  is  the  covariance  matrix  of  the  inputs  to  the  weights,  and  n  is  the 
number  of  weights.  Formula  (21)  applies  when  the  eigenvalues  are  all  equal. 
This  is  the  case  for  the  system  of  Fig.  2.  Multiplying  by  the  sampling  period  T, 
the  time  constant  is  expressed  in  seconds  of  real  time  as 


For  the  two-weight  adaptive  filter  of  Fig.  2, 


frace  R  =  ~CZ  +  —  Cz 
2  2 

=  C2 .  (23) 

This  is  the  sum  of  the  power  into  the  weights.  Combining  equations  (23),  (22), 
and  (19)  yields 


(BW)  =  ——  Hz. 

Tms« 


(24) 


Thus,  the  bandwidth  of  the  notch  is  the  reciprocal  of  the  time  constant  of  the 


learning  process,  for  the  simple  system  of  Fig.  2. 

Figure  5  shows  the  results  of  two  experiments  performed  to  demonstrate 
that  the  adaptive  system  acts  like  a  notch  filter.  In  the  first,  the  primary  input 
was  a  cosine  wave  of  unit  power  stepped  at  512  discrete  frequencies.  The  refer¬ 
ence  input  was  a  cosine  wave  with  frequency  u0  of  tt/2 T  rad/s.  The  value  of  C 
was  1.  and  the  value  of  fj.  was  1.25  x  10-2  .  The  frequency  resolution  of  the  fast 
Fourier  transform  was  512  bins.  The  output  power  at  each  frequency  is  shown 
in  Fig.  5(a).  As  the  primary  frequency  approaches  the  reference  frequency, 
significant  cancellation  occurs.  The  weights  do  not  converge  to  stable  values 
but  "tumble"  at  the  difference  frequency,*  and  the  adaptive  filter  behaves  like  a 
modulator,  converting  the  reference  frequency  into  the  primary  frequency. 
The  theoretical  notch  width  between  half-power  points.  1.59  x  10~2  u0.  compares 
closely  with  the  measured  notch  width  of  1.62  x  10~2  tJ„. 

In  the  second  experiment,  the  primary  input  was  composed  of  uncorrelated 
samples  of  white  noise  of  unit  power.  The  reference  input  and  the  processing 
parameters  were  the  same  as  in  the  first  experiment.  An  ensemble  average  of 
4096  power  spectra  at  the  noise  canceller  output  is  shown  in  Fig.  5(b).  An 
infinite  null  was  not  observed  in  this  experiment  because  of  the  finite  fre¬ 
quency  resolution  of  the  spectral  analysis  algorithm. 

In  these  experiments,  the  filtering  of  a  reference  cosine  wave  of  a  given 
frequency  caused  cancellation  of  primary  input  components  at  adjacent  fre¬ 
quencies.  This  result  indicates  that,  under  some  circumstances,  primary  input 

components  may  be  partially  cancelled  and  distorted  even  though  they  are  not 

«. 

correlated  with  the  reference  input.  In  practice  this  kind  of  cancellation  is  of 

concern  only  when  the  adaptive  process  is  rapid;  that  is,  when  it  is  effected  with 

•When  the  primary  and  reference  frequencies  are  held  at  a  constant  difference,  the  weights 
develop  a  sinusoidal  steady  state  at  the  difference  frequency.  In  other  words,  they  converge 
on  a  dynamic  rather  than  a  static  solution. 
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large  values  of  /z.  When  the  adaptive  process  is  slow,  the  weights  converge  to 
values  that  are  nearly  fixed,  close  to  the  Wiener  solution,  and  though  signal 
cancellation  as  described  in  this  section  occurs,  it  is  generally  not  significant 
due  to  the  fact  that  the  notch  is  extremely  narrow.  In  any  event,  the  primary 
input  appears  at  the  output  having  gone  through  a  notch  filter. 

V.  SIGNAL  CANCELLATION  PHENOMENA 

WITH  A  FROST  ADAPTIVE  BEAMFORMER 

Figure  6  shows  the  antenna  array  and  Frost  beamformer  that  were  used  in 
the  series  of  computer  simulation  experiments  to  be  presented  below.  The  sig¬ 
nal  was  assumed  to  be  incident  from  broadside  and  the  look  direction  con¬ 
straint  was  set  to  unit  gain  and  zero  phase  from  zero  frequency  to  half  the  sam¬ 
pling  rate,  i.e.  to  a  flat  response  over  all  frequencies.  The  jammer  was 
sinusoidal  at  one  quarter  the  sampling  frequency.  In  these  experiments, 
ambient  noise  and  receiver  noise  were  negligible.  A  typical  converged  beam 
pattern  is  shown  in  Fig.  6,  plotted  at  the  jammer  frequency.  Only  half  the  sym¬ 
metrical  pattern  is  plotted.  One  can  see  that  the  look  direction  gain  is  unity, 
and  that  the  gain  in  the  jammer  direction  is  .019,  34  dB  below  the  main  beam 
gain.  Observation  of  the  beam  pattern  gives  the  appearance  that  the  beam- 
former  is  working  perfectly. 

A  bandpass  signal  was  received  by  this  adaptive  beamformer  whose  spec¬ 
trum  is  shown  in  Pig.  7(a).  The  spectrum  of  the  sinusoidal  jammer  is  shown  in 
fig.  7(b).  The  output  signal  spectrum  is  shown  in  Fig.  7(c).  The  input  signal 
appears  at  the  output  having  gone  through  a  notch  filter.  The  notching  effect  is 
evident  in  the  output  signal  spectrum  and  is  indicative  of  gross  signal  distor¬ 
tion  at  the  beamformer  output.  The  notch  width  is  not  exactly  equal  to  the 
reciprocal  of  the  learning  curve  time  constant,  but  exceeds  it  by  a  factor  of  2. 
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The  conditions  for  the  derivation  of  the  notch  width  formula  (24).  i.e., 
sinusoidal  signals  appearing  with  exact  90°  separation  at  the  inputs  to  the  two 
weights,  are  not  precisely  met  with  the  16-weight  Frost  processor  under  the 
above  stated  experimental  conditions.  Nevertheless,  the  simple  formula  (24) 
does  give  at  least  an  approximate  prediction  of  notch  width  that  is  applicable  at 
most  jammer  angles. 

The  notching  phenomenon  in  adaptive  beamformers  is  somewhat  more 
complicated  than  in  adaptive  noise  cancelling  systems.  The  useful  signal  arriv¬ 
ing  from  the  look  direction  encounters  a  unit  gain  due  to  the  main  beam  con¬ 
straint.  This  is  analogous  to  the  direct  primary  signal  path  of  Fig.  1.  The  jam¬ 
mer  signal  arriving  at  other  than  the  look  direction  encounters  an  adaptive 
filter,  analogous  to  the  reference  signal  path  of  Fig.  1.  The  weights  in  the  adap¬ 
tive  beamf ormer  are  not  completely  free  as  they  are  in  the  adaptive  filter  of 
Fig.  1.  The  Frost  constraint  reduces  the  number  of  degrees  of  freedom  to  be 
equal  to  the  number  of  weights  multiplied  by  the  factor  ( [k~\)/k ,  where  k  is 
equal  to  the  number  of  antenna  elements.  Also,  the  spatial  processing  effects 
of  the  array  and  the  multichannel  structure  of  the  adaptive  processor  intro¬ 
duce  differences  in  the  dynamics  of  convergence  between  the  adaptive  beam- 
former  and  the  adaptive  noise  canceller. 

Additional  experiments  were  conducted  with  the  system  of  Fig.  6.  The  jam¬ 
mer  was  again  sinusoidal,  while  the  look-direction  signal  was  composed  of  white 
noise  of  unit  power.  The  jammer  power  was  varied.  Spectra  of  the  beamformer 
outputs  are  shown  in  Fig.  8.  With  the  jammer  power  set  at  its  lowest  level,  the 
signal  cancellation  notch  is  at  its  smallest  bandwidth  as  is  seen  in  Fig.  8a.  As 
the  jammer  power  is  increased  and  other  parameters  held  constant,  the  notch 
width  increases.  Figure  8(c)  shows  the  widest  notch  for  the  strongest  jammer 
signal  that  was  applied.  In  all  of  the  illustrated  cases,  the  relationship  between 
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notch  width  and  reciprocal  adaptive  time  constant  has  been  preserved. 

The  results  of  another  experiment  arc  shown  in  Pig.  9.  Here,  the  signal  was 
white,  and  the  jammer  was  a  strong  bandpass  noise.  Signal  components  were 
partially  cancelled  over  the  entire  jammer  spectral  band,  corresponding  to 
extensive  signal  distortion.  It  should  be  realized  that  results  of  this  type  would 
only  occur  in  cases  of  very  rapid  adaptation.  For  the  experiment  of  Fig.  9.  the 
time  constant  of  the  adaptive  process  was  approximately  equal  to  80  sampling 
periods.  The  bandwidth  of  the  jammer  was  approximately  equal  to  15%  of  the 
sampling  rate. 

Two  remedies  to  the  signal  notching  or  cancellation  problem  will  be  dis¬ 
cussed.  The  first  is  a  method  devised  by  K.  Duvall  based  on  the  use  of  two  signal 
processing  systems,  one  to  perform  the  adaptation,  the  other  to  generate  the 
system  output  signal. 


VI.  THE  DUVALL  BEAMFORMER 

In  order  to  prevent  signal  cancellation,  the  useful  signal  arriving  from  the 
look  direction  is  excluded  from  the  beamformer  in  which  the  adaptive  process 
takes  place.  In  the  system  shown  in  Fig.  10,  the  adaptive  process  is  only  used  to 
derive  a  set  of  weights.  These  weights  are  copied  into  a  separate,  identical 
"slaved"  processor  used  to  form  the  output  signal.  The  adaptive  process  could 
be  the  Frost  algorithm,  as  indicated  in  Fig.  10,  or  it  could  be  any  one  of  a 
number  of  adaptive  beamformer  algorithms  that  have  appeared  in  the  litera¬ 
ture. 

The  antenna  array  elements  in  Fig  10  are  assumed  to  be  uniformly  spaced 
along  a  line.  It  is  clear  from  the  block  diagram  shown  in  this  figure  that, 
because  of  the  subtractive  preprocessing,  the  look  direction  signal  will  not 
appear  at  the  Frost  adaptive  beamformer  inputs  but  that  the  jammer  signals 
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will  be  present  as  indicated.  Receiver  noise  and  ambient  noise  arc  neglected  in 
this  discussion.  The  Frost  algorithm  will  null  jammer  J.  The  strength  and  char¬ 
acter  of  signal  S  will  have  no  effect  on  the  weights.  Copying  the  weights  will 
cause  the  slaved  processor  to  have  a  main  beam  which  conforms  to  the  Frost 
constraints  established  for  the  look  direction,  and  to  have  a  null  in  the  exact 
direction  of  the  jammer  J.  The  correct  alignment  of  the  null  is  assured  since 
the  relative  phases  of  the  jammer  components  are  the  same  in  the  slaved 
beamformer  as  in  the  Frost  adaptive  beamformer,  where  the  nulls  originate. 
The  phasor  diagram  in  Fig.  10  verifies  this  in  accord  with  the  following  argu¬ 
ment. 

The  jammer  components  received  by  the  antenna  elements  are  indicated 
by  a  set  of  equal  amplitude  uniformly-spaced  phasors  Jo.  J 1.  J  z.  /a.  and  J+.  The 
phasor  inputs  to  the  Frost  beamformer  are  J \-Jo.  Jz~^v  Ja~J 2.  and  J4— J3. 
They  too  are  uniform-amplitude,  equally-spaced,  and  separated  by  the  same 
angles  as  the  received  jammer  components  J  q,  J Jz.  J 3.  and  J 4.  Since  the 
relative  phase  angles  are  the  same  in  the  slaved  processor  as  in  the  Frost  pro¬ 
cessor,  the  beam  pattern  notch  is  formed  at  the  proper  bearing  angle. 

The  Duvall  beamformer  uses  a  standard  adaptive  beamformer  as  one  of  its 
components.  Beam  steering  in  any  direction  can  be  accomplished  simply  by 
including  beam  steering  delays  in  the  antenna  circuits.  Phase  shifters  would 
be  adequate  for  narrowband  processes;  delay  lines  would  be  required  for  broad¬ 
band  processes.  It  should  be  noted  that  although  the  phasor  argument  applies 
only  to  one  jammer  at  one  frequency,  linearity  and  superposition  show  that  the 
principle  is  applicable  to  multiple  jammers  and  to  broadband  as  well  as  to  nar¬ 
rowband  jammers  and  signals. 

Experiments  have  been  performed  with  this  system,  and  results  are  given 
in  Figs.  11  and  12.  Figure  11  compares  the  output  spectrum  of  the  Frost  beam- 
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former  with  that  of  the  Duvall  beamformer  (using  a  Frost  beamformer),  both 
adapting  with  a  time  constant  of  .05,  with  the  same  array  and  with  the  same  sig¬ 
nal  and  jammer.  The  array  and  jammer  were  as  shown  in  Fig.  6.  After  perform¬ 
ing  the  comparative  experiments,  the  Frost  beamformer  showed  evidence  of 
strong  signal  cancellation,  while  the  Duvall  beamformer  showed  no  evidence  of 
signal  cancellation.  In  the  time  domain,  Fig.  12  compares  the  look-direction 
input  signal  with  the  output  signals  of  the  Frost  and  Duvall  beamformers.  In 
both  cases  the  weights  were  initialized  to  zero,  and  adaptive  transients  are  visi¬ 
ble  at  the  beginnings  of  the  output  tracings.  Beyond  the  region  where  the  tran¬ 
sient  exists,  substantial  signal  distortion  in  the  Frost  beamformer  output  is 
present.  The  distortion  power  was  measured  to  be  6  dB  below  the  input  signal 
power.  Such  distortion  is  not  apparent  at  the  output  of  the  Duvall  beamformer. 
Here  the  distortion  was  measured  to  be  110  dB  below  signal  level. 

The  Duvall  beamformer  appears  to  be  an  important  development  toward 
mitigating  the  effects  of  signal  cancellation.  It  is,  however,  a  recent  develop¬ 
ment,  and  possible  limitations  on  its  performance  have  yet  to  be  assessed. 
Effects  of  component  inaccuracies  and  array  imperfections  are  not  yet  under¬ 
stood.  How  to  use  it  with  other  than  straight  line,  evenly  spaced  arrays  is  not 
yet  determined.  Other  methods  for  eliminating  or  reducing  signal  cancellation 
effects  are  also  being  pursued,  such  as  spatial  dither  algorithms. 

VII.  SPATIAL  DITHER  ALGORITHMS 

Spatial  dither  algorithms  have  been  newly  conceived  for  the  purpose  of 
applying  locally  controlled  modulation  to  signals  arriving  at  angles  other  than 
the  look  direction  while  leaving  inputs  from  the  look  direction  unmodulated 
and  undistorted.  The  effect  is  to  cause  jammer  power  to  be  spread  spectrally, 
thereby  reducing  jammer  power  density.  When  used  with  a  conventional  adap¬ 
tive  beamformer,  spatial  dither  reduces  signal  cancellation  effects.  The  same 
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process  has  the  additional  capability  of  modulating  a  "smart"  jammer  signal  in 
a  way  that  is  totally  unpredictable  to  the  jammer,  thus  in  many  cases  render¬ 
ing  it  "less  smart." 

A  conceptually  simple  form  of  spatial  dither  algorithm  is  the  "3/4-inch  ply¬ 
wood"  approach,  pictured  in  Fig.  13.  The  elements  of  an  antenna  array  may  be 
imagined  to  be  fixed  to  a  piece  of  plywood  which  provides  a  rigid  insulating  sup¬ 
port,  so  that  the  entire  array  may  be  moved  mechanically.  In  either  one  or  two 
dimensions,  the  array  is  moved  in  directions  which  are  orthogonal  to  the  look 
direction.  Far-field  emanations  arriving  from  the  look  direction  will  be  undis¬ 
torted  by  the  mechanical  motion,  while  emissions  from  off  axis  sources  will  be 
distorted  by  an  unusual  shift-of-tim e-base  form  of  modulation.  (Electronic 
means  of  implementation  of  this  spatial  dither  process  are  being  devised.) 

The  outputs  of  the  antenna  elements  of  Fig.  13  could  be  applied  to  a  time 
delay  and  sum  (nonadaptive)  beamformer,  to  a  conventional  adaptive  beam- 
former.  or  to  a  Duvall  adaptive  beamformer.  Spatial  dither  could  be  beneficial 
in  each  case.  By  reducing  jammer  power  density,  some  antijam  protection  is 
provided  without  adaptive  beamforming,  and  additional  antijam  protection  is 
provided  with  adaptive  beamforming.  Reduction  of  signal  cancellation  effect  in 
a  Frost  beamformer  can  be  obtained  by  using  spatial  dither  preprocessing. 
Breakup  of  jammer  signal  structure  is  a  possible  form  of  signal  preprocessing 
applicable  to  all  types  of  adaptive  and  nonadaptive  beamformers. 

The  3/4-inch  plywood  approach  has  been  computer  simulated,  and  results 
are  presented  in  Fig.  14.  The  motion  was  random  and  was  done  along  a  line  per- 
pendicular  to  the  look  direction.  At  every  fourth  sample  time,  the  plywood  posi¬ 
tion  was  switched;  the  new  position  was  drawn  randomly  from  a  uniform  distri¬ 
bution  which  ranged  from  zero  to  eight  wavelengths.  Figure  14(a)  shows  the 
spectrum  of  the  look-direction  input  signal,  in  this  case  a  random  bandpass  sig- 


nal.  The  sinusoidal  jammer  spectrum  is  shown  in  Fig.  14(b).  The  spectrum  of 
the  jammer  from  the  physical  reference  frame  of  the  array  is  shown  in  Fig. 
14(c).  It  is  clear  that  the  jammer  power  is  greatly  spread,  that  jammer  power 
density  is  significantly  reduced,  and  that  the  jammer  signal  is  severely  dis¬ 
torted  from  its  original  form.  In  the  simulation,  bandpass  filters  were  used  with 
each  antenna  output  to  represent  the  effects  of  a  receiver  for  each  antenna. 
The  filtered  signals  w'ere  then  applied  to  a  conventional  Frost  adaptive  beam- 
former.  Some  signal  distortion  is  evident,  but  the  amount  of  distortion  is 
greatly  reduced  by  the  spatial  dither.  The  output  spectrum  shown  in  fig.  14(d) 
is  far  less  distorted  than  that  of  Fig.  7(c),  a  comparable  spectrum  obtained 
without  spatial  dither. 

VIII.  CONCLUSION 

Signal  cancellation  effects  occur  in  conventional  adaptive  beamformers 
when  jammer  power  and  adaptation  rate  are  high.  These  effects  can  cause  sig¬ 
nal  loss  in  the  case  of  narrowband  signals  or  cause  significant  signal  distortion 
in  the  case  of  wideband  signals.  Means  of  combatting  signal  cancellation  have 
been  proposed,  namely  the  Duvall  beamformer  and  the  spatial  dither  algorithm. 
The  latter  approach  will  probably  not  be  as  effective  as  the  former  against  sig¬ 
nal  cancellation  but  has  the  capability  of  destroying  "smart"  jammer  signals. 
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Results  of  two-weight  adaptive  noise  cancelling  experiments,  (a)  Pri¬ 
mary  input  composed  of  cosine  wave  at  512  discrete  frequencies,  (b) 
Primary  input  composed  of  uncorhelnted  samples  of  white  noise. 
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Power  spectra  of  Frost  beamformer  simulation,  (a)  Input  signal  spec¬ 
trum.  (b)  Jammer  spectrum,  (c)  Beamformer  output  spectrum. 
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Fxg.  B.  Frost  beamformer  output  spectra  for  white-noise  signal  of  power  1.0 
and  sinusoidal  jammer,  (a)  Jammer  power  =  12.5  (b)  Jammer  power  -• 
25.0  (c)  Jammer  power  =  50.0 
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Look-Direction  Signal.  S  /  Jammer. 


Fig.  10.  Duvall  beamformer  for  eliminating  cancellation  jamming. 


Fig.  11.  Comparison  of  Frost  and  Duvall  bearnformer.  (a)  Input  signal  spec¬ 


trum.  (b)  Jammer  spectrum,  (c)  Duvall  bearnformer  output  spec¬ 
trum.  (d)  Frost  bearnformer  output  spectrum. 


Fig.  12.  Comparison  of  time-domain  outputs  of  Frost  and  Duvall  beamformers. 

(a)  Signal  input,  (b)  Duvall  bearnformer  output,  (c)  Frost  bearnformer 
'  output. 
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ADAPTIVE  POWER  SEPARATORS 


Timothy  Saxe 


1  --  INTRODUCTION 


( 


An  adaptive  power  separator  (APS)  is  an  adaptive  filter  that 
discriminates  between  signals  on  the  basis  of  their  powers.  There  are 
different  varieties  of  APS.  Some  pass  low  power  signals  and  reject  high 
power  signals  (low-power  pass  APS),  some  pass  only  high  power  signals  * 
(high-power  pass  APS),  and  others  pass  only  signals  with  power  within  a 
specified  range.  A  low-power  pass  APS  could  be  used  in  the  "front  end" 
of  a  radio  receiver  to  eliminate  powerful  signals  from  nearby  trans¬ 
mitters  while  still  passing  weak  signals  from  distant  transmitters.  Of 
course,  a  single  channel  APS  cannot  distinguish  between  two  signals  that 
overlap  spectrally. 

This  paper  develops  and  analyses  three  different  single  channel 
low-power  pass  adaptive  power  separators.  Therefore,  in  this  paper  the 
term  adaptive  power  separator  is  understood  to  mean  a  single  channel 
power  separator  that  passes  only  low  power  signals.  However,  the  power 
separators  developed  in  this  paper  can  be  easily  altered  to  pass  high 
power  signals  or  even  signals  in  a  specified  power  range. 

Previous  workers  in  the  field  of  adaptive  filtering  [Wid,  Zah, 

Tre,  Che]  have  noted  that  some  adaptive  filters  have  a  property 
called  power  inversion.  Power  inversion  is  a  weak  form  of  power  separa¬ 
tion  in  which  signals  with  high  power  at  the  input  to  the  adaptive  system 
Tiave  less  output  power  than  signals  with  low  input  power.  However,  the 
transition  from  low  gain  to  high  gain  is  not  sharp  so  these  filters 
cannot  discriminate  between  two  signals  that  have  nearly  the  same  power. 
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The  purpose  of  this  work  is  to  develop  an  APS  with  a  sharp  transition  ^ 

that  will  be  able  to  discriminate  between  two  signals  with  nearly  the 
same  power. 

Any  APS  is  a  nonlinear  system.  However,  all  of  the  adaptive  power 
separators  described  in  this  work  become  linear  filters  if  the  adaptive 
process  is  stopped.  The  adaptive  process  produces  linear  filters  that 
use  frequency  discrimination  to  reject  high  power  signals.  Consequently, 
a  single  channel  APS  cannot  discriminate  between  two  signals  that  over¬ 
lap  in  frequency.  This  limitation  may  be  overcome  by  using  a  multi¬ 
channel  APS  if  the  signals  are  spatially  separated.  For  example,  an  APS 
built  into  an  antenna  array  could  discriminate  between  two  signals 
provided  that  they  either  had  different  angles  of  arrival  or  they  did 
not  overlap  in  frequency. 

c 


C 
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2  --  PREVIOUS  WORK:  THE  LALE  APS 

The  single  channel  adaptive  system  which  is  called  the  adaptive 
(spectral)  line  enhancer  (ALE)  was  first  proposed  by  Widrow  [Wid2] 
for  the  detection  of  sinusoidal  components  buried  in  noise.  He  showed 
that  the  gain  of  the  ALE  was  a  function  of  the  signal -to-noi se  ratio. 
Thus  the  ALE  can  be  used  as  an  APS.  However,  the  threshold  power*  of 
the  ALE  is  difficult  to  control  since  it  is  a  function  of  the  input 
noise  power  as  well  as  the  signal  power.  Widrow  and  Treichler  [Wid3, 
Wid4,  Tre]  proposed  a  modified  ALE,  called  the  ''leaky"  ALE  (LALE),  that 
algorithmically  simulates  the  effect  of  added  input  noise  without 
actually  adding  noise  to  the  input.  By  this  means  the  effective  signal- 
to-noise  ratio,  and  hence  the  threshold  power  of  the  LALE,  can  be 
controlled  by  altering  a  parameter  of  the  adaptive  process.  Thus  the 
LALE  is  an  APS  with  a  selectable  threshold  power. 

In  this  chapter  we  will  analyze  the  performance  of  the  LALE  and 
show  that  it  can  be  used  as  a  either  a  high-power  pass  APS  or  a  Tow- 
power  pass  APS.  A  brief  intuitive  explanation  of  the  LALE  is  presented 
first  to  aid  in  understanding  the  mathematical  analysis  of  the  LALE. 

As  was  mentioned  earlier,  the  ALE  can  be  used  as  an  APS.  Since 
the  ALE  is  simpler  to  understand  than  the  LALE,  we  will  first  study  the 
operation  of  the  ALE.  Figure  2-1  shows  a  block  diagram  of  an  ALE. 

The  heart  of  the  ALE  is  the  LMS  (Least  Mean  Square)  adaptive  filter 

*The  threshold  power  of  an  APS  is  defined  as  the  signal  power  for 
which  the  voltage  gain  of  the  APS  is  1/2. 
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[Wid].  An  LflS  filter  is  a  finite  impulse  response  filter  that  tries  to 
form  the  best  least  squares  fit  between  its  input  signal  (x(k)  in  the 
ALE)  and  the  primary  input  (d(k)  in  the  ALE).  The  "error’1  output, 
which  is  the  difference  between  the  primary  input  and  the  adaptive 
filter’s  output,  measures  how  successful  the  filter  has  been  in  matching 
its  input  to  the  primary  input.  The  "error"  output  is  used  as  an 
auxiliary  input  to  the  LMS  filter  so  that  the  adaptive  algorithm  can 
decide  how  to  modify  the  transfer  function  of  the  filter  to  reduce  the 
mean  square  error. 

We  can  analyze  the  performance  of  the  ALE  without  knowing  how  the 
LMS  filter  actually  reduces  the  mean  square  error.  For  the  sake  of 
simplicity,  assume  that  the  input  to  the  ALE  consists  of  a  sinusoid  in 
white  noise.  The  input  to  the  adaptive  filter,  x(k),  is  a  delayed 
version  of  the  primary  input,  d(k).  Therefore  the  "error"  signal,  e(k), 
is  the  sum  of  three  terms:  1)  noise  in  the  primary  input,  2)  noise 
passing  through  the  adaptive  filter,  and  3)  a  sinusoidal  component.  The 
significance  of  the  delay  between  d(k)  and  x(k)  is  that  the  adaptive 
filter  cannot  match  the  noise  in  its  input  with  the  noise  in  the  primary 
input,  but  can  match  the  sinusoidal  component  by  generating  the  appro¬ 
priate  gain  and  phase  at  the  frequency  of  the  sinusoid.  To  minimize 
the  mean  square  error,  the  adaptive  filter  will  try  to  pass  as  little 
noise  as  possible  while  also  passing  the  sinusoid  with  the  appropriate 
gain  and  phase. 

A  filter  with  a  sinusoidal  weight  vector,  at  the  same  frequency 
as  the  input  sinusoid,  will  have  a  bandpass  characteristic  at  the 
frequency  of  the  input  sinusoid  and  minimum  bandwidth.  This  is  exactly 
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low  SNR,  i.e.  low  signal  power  for  fixed  noise  power,  the  gain 
a)  is  low,  but  for  high  SNR  the  gain  tends  to  one.  Thus  the  ALE 
is  a  high-power  pass  APS  because  it  passes  strong  signals  better  than 
weak  signals.  Figure  2-2  shows  the  transfer  function  of  an  ALE  for 
several  different  SNR's.  Incidental ly,  Treichler  [Tre]  has  shown  that 
the  response  to  several  sinusoidal  inputs  is  the  superposition  of  the 
response  to  the  individual  sinusoids  provided  that  the  frequencies  of 
all  the  sinusoids  differ  by  more  than  1/nT  Hz  (where  T  is  the  sampling 
period). 

A  deficiency  of  the  ALE  is  that  the  response  depends  upon  the 
signal  to  noise  ratio,  not  the  signal  power.  Widrow  and  Treichler 
[Wid3]  proposed  adding  noise  to  the  input  of  the  ALE,  which  would  allow 
the  designer  to  control  the  signal  to  noise  ratio  and  thereby  control 
the  characteristics  of  the  ALE.  Of  course  this  technique  makes  the  out¬ 
put  of  the  ALE  more  noisy  and  increases  the  variance  of  the  LMS  filter's 
weights,  which  modulates  the  output  signal.  However,  since  the 
characteristics  of  the  injected  noise  are  known,  the  effects  of  the 
noise  upon  the  LMS  algorithm  can  be  computed.  These  effect  can  then  be 
incorporated  into  the  LMS  algorithm  [Gri,  Wid4,  Tre]  which  results,  in 
the  case  of  white  noise  added  to  the  input,  in  the  leaky  LMS  (LLMS) 
algorithm.  The  advantages  of  using  the  LLMS  algorithm  instead  of 
injecting  noise  are  a  less  noisy  output,  less  weight  variance,  no  need 
to  build  a  noise  generator  and  the  LLMS  algorithm  can  act  as  if  noise 
with  negative  power  were  added  to  the  input.  For  convenience  we  will 
refer  to  "algorithmic  noise"  when  we  mean  the  noise  effects  that  are 
incorporated  into  the  LLMS  algorithm. 
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Figure  2-2.  Transfer  function  of  the  ALE 
for  three  powers  of  an  input 
signal 
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By  altering  a  parameter  of  the  LLMS  filter  the  designer  can  control 
the  power  of  the  "algorithmic  noise."  This  permits  the  designer  to 
control  the  effective  SNR  (effective  SNR  is  the  signal  power  divided  by 
the  sum  of  the  actual  input  noise  power  and  the  "algorithmic  noise" 
power),  and  thereby  control  the  characteristics  of  the  LALE.  Thus  the 
LALE  is  a  high-power  pass  APS  with  a  controllable  characteristic. 

The  rest  of  this  section  is  devoted  to  a  more  comprehensive  and 
quantitative  analysis  of  the  LALE.  This  analysis  is  based  on  Treiclfler's 
[Tre]  extensive  analysis  of  the  ALE  with  appropriate  modifications  to 
account  for  the  changes  between  the  LMS  (Least  Mean  Square  [Wid]) 
algorithm  used  in  the  ALE  and  the  LLMS  (Leaky  LMS)  algorithm  used  in 
the  LALE.  While  the  analysis  is  for  the  most  part  general,  assumptions 
about  the  inputs  are  required  to  obtain  insight  into  the  solutions  for 
the  LALE  weight  vector.  In  particular,  the  LALE  APS  can  be  used  to 
reject  high-power,  narrowband  interference  (where  narrowband  means  a 
bandwidth  significantly  less  than  4jr  where  n  is  the  number  of  LALE 
weights  and  T  is  the  sampling  interval).  Therefore  we  assume  that  the 
input  to  the  LALE  consists  of  narrowband  signals  and  white  noise. 
Treichler  [Tre]  showed  that  under  these  conditions  the  input  to  the 
LALE  can  be  modeled  as  one  sinusoid  per  narrowband  signal.  Thus  the 
input  to  the  LALE  is  modeled  as  a  number  of  sinusoids  in  white  noise. 
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We  start  by  defining  various  quantities  in  accordance  with  figures 


( 


2-3  and  2-4; 

x(k)  =  adaptive  filter  input  at  time  k 
X(k)  =  [x(k),  x(k-l),  ....  x( k-n+1 )]T 

=  contents  of  tapped  delay  line  at  time  k 
W(k)  =  vector  of  n  filter  weighting  coefficients 
y(k)  =  XT(k)W(k) 

'  =  adaptive  filter  output 
e(k)  =  d(k)-y(k) 

=  the  "error"  signal  (drives  the  adaption) 
d(k)  =  x(k+A) 

=  desired  (or  reference)  signal 
&  =  decorrelation  delay  time 

The  LLMS  algorithm  for  modifying  the  weight  vector  is: 

W( k+1 )  =  vW(k)  +  2ye(k)X(k) 

=  [vl  -  2yX( k)XT( k)]W(k)  +  2yd(k)X(k)  (2-1) 

where 

v  =  an  arbitrary  constant  (leak  factor) 

=  1  -  2  py 

Y  =  power  of  "algorithmic  noise” 

y  =  arbitrary  constant  which  controls  the  rate  of  adaption  and 
the  stability  of  the  algorithm 

( 
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Figure  2-4.  Diagram  of  the  Leaky  Adaptive  Line  Enhancer. 


The  average  behavior  of  the  LALE  weight  vector  can  be  derived  from 
equation  2-1  by  taking  the  expected  value  of  both  sides  of  Eq.  2-1  and 
ignoring  the  covariance  between  X(k)XT( k)  and  W(k)  (this  is  an  acceptable 
approximation  if  p  is  small  [Sen,  Dan]); 

E{H(k+l)}  *  [vl  -  2yR(k)]  E{U(k)}  +  2yP(k)  (2-2) 

where 

=  E{x(k-i)  x(k-j)} 

=  input  autocorrelation  matrix. 

CP(k)3i  *  E{d(k)  x(k-i)} 

=  desired-to-input  cross 
correlation  vector 

If  the  input,  x(k),  consists  of  one  or  more  statistically 

independent,  zero-mean  and  stationary  components,  then  R(k)  =  Rn+Rs 

and  P{k)  =  p  +P  ,  where  R  is  the  autocorrelation  matrix  associated 
'  '  n  s  n 

with  the  inputs  whose  correlation  times  are  less  than  a  and  Rs  is  the 

autocorrelation  matrix  associated  with  the  inputs  whose  correlation 

times  are  greater  than  a.  If  we  assume  that  the  incoherent  components 

2  2 

can  be  modeled  as  white  noise  with  power  a  ,  then  Rn  =  a  I.  Finally, 
if  A>.  1  then  P(k)  =  P$  since  Pn  is  equal  to  zero  (it  is  the  autocor¬ 
relation  of  white  noise  for  lags  greater  than  A).  Provided  that  all  of 
these  conditions  hold,  Eq.  2-2  can  be  written  as: 

E{W( k+1 ) }  -  [vl  -  2u(a2I  +  Rs)]  E{W(k) }  +  2UPS  .  (2-3) 

Since  R$  is  the  autocorrelation  matrix  of  a  real  valued  stochastic 
process,  it  is  symmetric,  Toeplitz,  positive  semi-definite  and  has  a 
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full  set  of  n  orthonormal  eigenvectors  Qq  through  Qn  -j  [Gra].  These 
eigenvectors  can  be  used  to  form  the  modal  matrix  Q,  where 
Q  =  [Q0.  •••»  Qn_-| 3 •  It  is  well  known  [Der]  that  the  modal  matrix  has 
the  following  properties: 

QTR$Q  =  A  (2-4a) 

QTQ  =  I  (2-4  b) 

where  A  is  the  diagonal  matrix  of  the  eigenvalues  of  R$ 

c  denotes  the  complex  conjugate  of  c 

If  a  transformed  weight  vector  is  defined  as  W'(k)  =  Q^W(k)  and  if  both 
sides  of  Eq.  2-3  are  premultiplied  by  QT,  then: 

E{lT(k)}  =  QTE{W(k)} 

-  QT[vI  -  2p(o2I+Rs)]  QE{VT(k) }  +  2yQTPs  .  (2-5) 

Equation  2-5  can  be  further  simplified  by  using  the  relations  given  in 
Eq.  2-4  and  by  defining  C  =  QT  Ps: 

E(W'(k)}  =  [vl  -  2y(o2I  +  A)]  E{W'(k) }  +  2yC  .  (2-6) 

The  transition  matrix  of  this  recursion  expression  is  diagonal 
since  I  and  A  are  diagonal.  Therefore  the  transformed  weights,  W"(k), 
are  uncoupled,  so  a  separate  scalar  recursion  expression  can  be  written 
'•for  each  [E{W'(k)}]^.  The  expression  for  the  ith  uncoupled  weight  is: 

[E(W'(k)}]i  -=  [v  -  2y(o2+Xi)][E{W'(k)}]i  +  2y[C]i 

where  A .j  is  the  ith  eigenvalue  of  R$. 


(2-7) 


Incidentally,  if  v  =  1  -  2yy  is  substituted  into  Eq.  2-6,  the 
result  is: 

E{W'(k)}  =  [I  -  2y([a2+y]I  +A)]  E{ir(k}}  +  2yC  .  (2-8) 

Equation  2-8  shows  that  y  affects  the  average  weight  vector  in  the 
same  my  as  a  .  Since  the  effects  of  y  are  analagous  to  the  effects 
of  the  white  noise  in  the  input,  y  is  called  the  "algorithmic  noise"  . 
power.  An  LLMS  filter  with  parameter  y  will  have  the  same  average 
weight  vector  as  a  Wiener  filter  with  noise  of  power  y  added  to  its 
input  signal.  In  a  similar  vein,  the  term  "effective -noise"  power  is 
used  to  mean  [a  +  y],  which  is  the  sum  of  the  actual  input  noise  power 
and  the  "algorithmic  noise"  power. 

Using  the  definition  of  y,  the  scalar  recursion  for  the  uncoupled 
weights  (Eq.  2-7)  can  be  written  as: 

[E{VT(k)}].j  =  [1  -  2y(a2n+Xi)][E{W-(k)}]i  +  2u[C].  .  (2-9) 

If  the  adaption  constant  y  is  chosen  so  that 

0  <  |1  -  2y(a2+y+Xi) |  <  1  (2-10) 

then  each  uncoupled  recursion  expression  (Eq.  2-9)  is  stable,  and 
E(VT(k)}  will  converge  to  a  value  E(W'(«)}  which  is  independent  of  the 
starting  conditions.  The  converged  value  is: 

[CL 

[E{vr(~)}L  =  - 1 -  i=0,  1 . n-l  .  (2-11) 

o  +  y  +  Xi 

Reversing  the  coordinate  transformation  yields: 

n-]  [CD  - 

£(«(»)}  *  Z  -? - 1 - Q<  •  (2-12) 

i=0  o  +  y  +  Xj 
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Equation  2-12  defines  the  expected  value  of  the  converged  weight 
vector.  This  definition  is  in  terms  of  the  eigenvectors  of  the  input 
autocorrelation  matrix,  the  power  of  the  true  input  noise,  the  power  of 
the  equivalent  noise  and  the  power  of  the  coherent  signal  component. 

To  determine  the  response  of  the  LALE  to  particular  inputs  requires 
explicit  solutions  for  the  and  [C]^,  and  knowledge  of  a  and  y* 

Since  the  input  to  the  LALE  can  be  modeled  as  multiple  sinusoids  in 
white  noise,  Treichler's  results  [Tre]  are  useful.  Treichler  showed 
that  if  the  input  consists  of  one  or  more  sinusoids  which  meet  the 
following  criteria:  1)  frequencies  are  between  20%  and  80%  of  the 
Myquist  frequency  (otherwise  the  two  active  eigenvalues  degrade  into  a 
single  active  eigenvalue),  and  2)  all  sinusoids  are  separated  in 
frequency  by  at  least  Hz  (otherwise  they  interact  instead  of  acting 
independently),  then  for  each  sinusoid  of  power  p^  and  frequency  the 
following  approximations  hold: 


*1 
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These  results  can  be  substituted  into  Eq.  2-12  to  determine  the  average 

weight  vector  of  the  LALE.  Specifically,  if  the  input  to  the  LALE 

consists  of  a  single  sinusoid  of  power  p  and  frequency  u  in  white  noise 
2 

of  power  a  ,  then 

[E{W(«)}] -  =  -2— £ - pp  Cos[(i+A)uT]  ..  (2-14) 

°  +  y  +  2? 

A  filter  with  this  weight  vector  has  a  bandpass  characteristic  with  a 

n 

2? 

maximum  gain  of  — ? — — - — —  at  the  frequency  of  the  sinusoid-  Figure 
o  +  Y  +  JP 

2-5  shows  the  transfer  function  of  the  LALE  for  several  input  powers. 

The  operation  of  the  LALE  is  more  conveniently  described  as  a  gain 
characteristic  (see  Figure  2-6).  The  gain  characteristic  shows* the 
converged  gain  of  the  adaptive  system  at  the  frequency  of  the  input 
signal  as  a  function  of  the  signal's  power.  We  see  from  Figure  2-4  that 
the  LALE  passes  powerful  signals  and  rejects  weak  signals  and  is  thus  a 
high-power  pass  APS.  If  a  low-power  pass  APS  is  required,  the  "error" 
output  of  the  LALE  can  be  used  instead  of  the  "y"  output.  The  "error" 
output  gain  characteristic  is  shown  in  Figure  2-7. 

The  threshold  power  of  an  APS  hereby  defined  as  the  input  signal 

power  for  which  the  gain  of  the  APS  is  j.  In  the  case  of  the  LALE  the 

2  2  2 
threshold  power  is  -{cr  +y).  Now  the  input  noise  power  a  ,  is  beyond 

*  the  designer's  control,  but  y  is  controllable  since  it  depends  on  y  and 

v  which  are  both  parameters  of  the  LLMS  filter.  Thus  the  threshold 

power  of  the  LALE  is  easily  controllable.  However,  the  transition  from 

high  gain  to  low  gain  is  gradual  so  the  LALE  is  said  to  have  a  soft 
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Figure  2-5.  Transfer  function  of  the  LALE 
for  three  powers  of  an  input 
signal . 
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Figure  2-6.  "y"-gain  characteristic  of  LALE 
(High-power  pass  APS). 
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Figure  2-7.  "error"-gain  characteristic  of  LALE 
(Low-power  pass  APS). 
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threshold  or  broad  transition  region.  Because  of  the  soft  threshold  a 
LALE  based  APS  cannot  decisively  discriminate  between  two  signals  with 
nearly  the  same  power. 


( 

3  --  PRACTICAL  SPECIFICATIONS  FOR  AM  ADAPTIVE  POWER  SEPARATOR. 

Figure  3-1  shows  the  gain  characteristic  of  leaky  adaptive  line 
enhancer  superimposed  on  the  gain  charac ter istic  of  an  ideal  APS.  The 
gain  characteristic  of  the  ideal  APS  has  unit  gain  for  input  power  less 
than  some  threshold  power  and  zero  gain  for  input  power  greater  than 
the  same  threshold.  In  contrast,  the  gain  characteristic  of  the  LALE 
(leaky  adaptive  line  enhancer)  has  no  sharp  discontinuity  and  is  never 
zero.  Thus  the  LALE  cannot  reject  signals,  it  can  only  attenuate  them. 
However,  the  LALE  has  a  gain  of  nearly  one  for  low  power  inputs  and  a 
gain  of  nearly  zero  for  high  power  inputs.  Thus  the  LALE  is  similar 
to  the  ideal  APS  because  it  passes  low  power  signals  more  strongly 
than  it  passes  high  power  signals.  The  purpose  of  this  chapter  is  to 
define  specifications  which  quantify  the  difference  between  a  practical 
APS  (such  as  the  LALE)  and  the  ideal  APS. 

One  approach  to  defining  a  practical  gain  specification  is  to  alloy; 
a  range  of  gains  v;hich  are  near  to  the  ideal  gain.  For  example,  instead 
of  requiring  that  the  gain  be  one  if  the  input  power  is  less  than  some 
threshold,  require  that  the  gain  be  between  1-5^  and  l+d|  if  the  input 
power  is  less  than  the  threshold.  Figure  3-2  shows  a  gain  tolerance 
scheme  which  is  based  on  this  idea.  A  mathematical  description  of  the 
allowable  gains  is: 

1  -  61  <  1  +  51  for  P<P] 

g  <  1  +  for  p1  <  p  <  p2  (3-1) 

g  <  iSj  for  p2  <  p 
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We  will  refer  to  this  tolerance  schane  as  tolerance  scheme  one  (TSl). 

If  61  =  0,  62  =  0  and  P1  =  p2*  then  TSl  describes  the  ideal  APS  gain 
characteristic. 

An  alternative  is  to  consider  the  APS  entirely  from  the  power 
domain.  Figure  3-3  shows  a  power-out  vs.  power-in  curve  for  an  ideal 
APS.  The  power-out  vs.  power-in  formulation  suggests  bounding  of  the 
output  power  rather  than  the  filter  gain. 

A  practical  power-out  vs.  power-in  specification  can  be  derived 
by  permitting  tolerances  about  the  ideal  curve.  Figure  3-4  shows  a 
practical  power-out  vs.  power-in  tolerance  scheme.  This  tolerance 
scheme  will  be  called  tolerance  scheme  two  (TS2).  TSl  and  TS2  differ 
because  TS2  controls  the  response  to  high  power  inputs  more  tightly 
than  TSl.  If  TSl  is  converted  into  a  power-out  vs.  power-in  specifica¬ 
tion  (see  Figure  3-5),  it  is  clear  that  TSl  permits  high  power  inputs 
to  have  high  output  power.  TS2,  on  the  other  hand,  guarantees  that 
signals  with  input  power  greater  than  p2  will  have  less  output  power 
than  signals  with  an  input  power  between  Pg  and  p^. 

A  mathematical  description  of  TS2  is: 


1  -  Pi  <  9  <  1  +  P] 

for 

P  1  Pl 

g  <  1  +  Pl 

for 

Pl  < 

P  <  P2 

(3-2) 

fz 

9  i  p2  Jt 

for 

P2< 

P 

If  p1  =  0,  p2  =  0  and  p1  =  p2>  then  T$2  describes  the  ideal  APS 
characteristic.  Also,  if  =  pj  and  <5 £  =  P2  ^en  TS2  is  equivalent 
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to  TSl  for  input  power  less  than  p2  (but  TS2  is  more  stringent  than  TSl 
for  input  powers  greater  than  P2). 

The  two  tolerance  schemes  presented  in  this  chapter  can  be  used  to 
describe  a  practical  APS.  By  using  these  schenes  to  characterize 
different  power  separators  we  will  be  able  to  identify  the  strengths  and 
weaknesses  of  the  different  power  separators  with  the' goal  of  developing 
improved  practical  power  separators. 


4  --  APS2:  A  DPT  BASED  ADAPTIVE  POWER  SEPARATOR 


A  straightforward  approach  to  adaptive  power  separation  is  to 
evaluate  the  power  spectrum  of  the  input  signal  and  then  to  use  the 
spectrum  to  design  a  filter  which  has  zero  gain  at  frequencies  where  the 
power  spectral  density  exceeds  a  chosen  threshold  and  unit  gain  else¬ 
where.  When  the  input  signal  is  applied  to  the  resulting  filter 
(see  Figure  4-1),  the  output  will  contain  only  the  weak,  input  components 
which  do  not  spectrally  overlap  the  strong  input  components.  This  is  a 
general  structure  for  an  APS. 

The  power  spectrum  of  a  signal  can  be  evaluated  by  the  digital 
Fourier  transform  (DPT).  This  technique,  which  is  called  the  Welch- 
periodogram  method  [Opp] ,  creates  a  power  spectral  estimate  by  averaging 
the  squared  magnitude  of  the  DFT  of  the  input  data.  An  associated 
filtering  technique  multiplies  the  outputs  of  the  DFT  (which  correspond 
to  spectral  components)  by  zero  or  one  and  then  takes  an  inverse  DFT  to 
obtain  the  output  time  series.  Figure  4-2  shows  the  block  diagram  of  an 
adaptive  power  separator  which  is  based  on  these  ideas.  This  APS  is 
called  the  APS2. 

The  input  signal  is  applied  to  the  APS2  as  a  series  of  data  samples, 
and  the  APS2  computes  the  DFT  of  the  last  n  data  samples: 

n-1 

[G(k)]m  =  £  x(k-i)  Exp(27T^p)  m  =  0,  1 ,  ....  n-1  (4-1) 
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Figure  4-2,  Block  diagram  of  DFT/threshold  based  power  separator  (APS2) 


where 


G(k)  =  vector  of  OFT  values  at  time  k 
m  =  spectral  bin  number 
=  n  f  T 

T  *  sampling  period 
f  =  frequency 

Each  DFT  value  [G(k)3m  is  the  output  of  a  finite  impulse  response  (FIR) 
filter  with  transfer  function: 

n-1 

H(Z,m)  =  Z  Exp(j27T^-)Z'1 
i=0 

=  - k - m  =  0,  1,  ....  n-1  .  (4-2) 

1  -  Exp(j2^)Z_1 

Thus  the  DFT  of  the  input  sequence  can  be  computed  by  a  bank  of  n  FIR 
filters  with  transfer  functions  given  by  Eq.  4-2  [Pap].‘  Figure  4-3 
shows  an  implementation  of  the  DFT  which  is  based  on  Eq.  4-2.  Uhile 
a  more  efficient  way  to  compute  the  DFT  is  by  one  of  the  fast  algorithms 
such  as  the  fast  fourier  transform,  great  insight  can  be  obtained  by 
viewing  the  DFT  as  a  bank  of  bandpass  filters.  Thus  we  will  analyze  the 
APS2  as  if  it  were  made  from  a  bank  of  bandpass  filters. 

The  APS2  filtering  technique  multiplies  the  outputs  of  the  DFT  by 
zero  or  one  and  then  inverse  DFT’s  the  resulting  filtered  signal  to 
obtain  a  time  series.  If  the  DFT  is  modeled  as  a  bank  of  bandpass 
filters,  as  shown  in  Figure  4-3,  then  the  APS2  filtering  technique  can 
be  modeled  as  multiplying  the  outputs  of  the  bandpass  filters  by  zero  or 
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one  and  then  summing  the  weighted  outputs.  This  is  show n  in  Figure  4-4 
where  the  threshold  devices  have  an  output  of  1  if  the  input  power  at  a 
given  spectral  bin  is  less  than  the  threshold  power  and  0  if  the  input 
power  is  greater  than  the  threshold  power.  We  will  use  this  model  to 
analyze  the  two  phases  of  the  APS2's  operation:  spectral  estimation  and 
filter  design. 

Recall  that  the  APS2  uses  the  average  squared  magnitude  of  the  DFT 
as  an  estimate  of  the  power  spectrum.  Also,  the  DFT  output  is  the 
output  of  a  filter  with  transfer  function  given  by  Eq.  4-2.  This  filter 
has  n  zeros,  equally  spaced  on  the  unit  circle  at  frequencies  which  are 
integer  multiples  of  ^r,  and  a  pole  at  Exp(j2m2-)  which  cancels  the  zero 
at  frequency  The  transfer  function  of  such  a  filter  is  shown  in 

Figure  4-5  (m  =  4,  n  =  16  for  this  filter).'  If  the  input  consists  of 
sinusoids  at  integer  bin  numbers*  (i.e.  their  frequencies  are  integer 
multiples  of  --) ,  then  the  output  power  of  the  mth  filter  will  be  propor¬ 
tional  to  the  power  of  the  sinusoid  at  frequency  ^  because  the  zeros 
of  the  comb  filter  will  cancel  all  of  the  sinusoids  except  the  sinusoid 
at  frequency  where  there  is  a  pole  zero  cancellation.  Thus,  provided 
that  the  input  consists  of  on-bin  sinusoids,  the  power  in  the  OFT  out¬ 
put  is  the  true  power  spectrum.  However,  if  the  input  has  off-bin 
sinusoids,  then  the  APS2  will  not  estimate  the  true  power  spectrum. 
Consider  the  case  of  a  sinusoid  at  frequency  f  (see  Figure  4-5).  This 
sinusoid  is  in  a  sidelobe  of  every  DFT  filter  and  so  will  pass  through 

*In  digital  filtering,  signal  frequencies  are  often  expressed  in 
terms  of  spectral  bin  numbers  where: 

bin  number  =  f  x  T  x  number  of  points  in  the  DFT  . 

For  convenience,  signals  with  integer  bin  numbers  are  called  on-bin  . 
signals  and  those  with  non-integer  bin  numbers  are  called  off-bin 
signals. 
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In  summary,  the  APS2  functions  ideally  for  any  number  of  sinusoidal 
inputs  provided  that  their  frequencies  all  correspond  to  integer  bin 
numbers.  Sinusoidal  inputs  which  are  not  on-bin  can  cause  undesirable 
effects.  First,  the  DFT  based  power  spectrum  estimator  will  exhibit 
leakage.  This  causes  the  power  spectral  estimate  to  show  more  power 
than  is  really  present  at  seme  frequencies,  so  the  APS2  may  block  weak 
signals  because  the  power  spectral  estimate  falsely  indicates  that  they 
are  strong.  Similarly,  the  power  spectral  estimate  may  show  less  power 
at  some  frequencies  than  is  actually  present.  Second,  the  filter  design 
used  by  the  APS2  to  reject  powerful  signals  unfortunately  creates 
relatively  large  gains  in  portions  of  the  rejection  region.  Consequently 
the  APS2  will  only  partially  reject  powerful  off-bin  signals. 

Simulation  of  the  APS2 

An  APS2  based  on  a  32  point  DFT  was  simulated.  The  first  test 
input  to  the  APS2  was  the  sum  of  two  sinusoids.  Sinusoid  one  was  at  25% 
of  the  Nyquist  frequency  (bin  number  =  4)  and  sinusoid  two  was  at  56.25 % 
of  the  Nyquist  frequency  (bin  number  =  9).  During  the  simulation  the 
power  of  the  first  sinusoid  was  varied  from  0.1  to  100  and  the  power  of 
the  second  sinusoid  was  held  constant  at  0.1.  The  threshold  power  was  1. 

Figures  4-7A  and  B  show  the  gain  characteristic  of  the  APS2  at  the 
frequency  of  each  sinusoid.  Since  both  of  the  signals  were  on-bin  the 
DFT  did  not  suffer  from  leakage  and  the  APS2  functioned  perfectly. 

Figure  4-7A  shows  that  the  gain  at  the  frequency  of  sinusoid  one  was 
unity  v/hen  the  power  was  less  than  the  threshold  and  zero  when  the  power 
exceeded  the  threshold.  Figure  4-7B  shows  that  the  gain  at  the  frequency 
of  the  sinusoid  two  remained  unity  for  all  powers  of  sinusoid  one. 
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Figures  4-8  shows  the  transfer  function  of  the  APS2  for  several 
powers  of  signal  one.  from  this  figure  we  see  that  the  APS2  cancels 
the  powerful  signal  by  creating  a  notch  at  its  frequency. 

Figures  4-9A  and  3  show  the  gain  of  the  APS2  for  the  another  test 
input.  Again  the  test  input  was  the  sum  of  two  sinusoids.  Sinusoid  one 
v/as  at  28.125%  of  the  ilyquist  frequency  (bin  number  -  4.5)  and  sinusoid 
two  was  at  56.25%  of  the  Nyquist  frequency  (bin  number  -  9).  During  the 
simulation  the  power  of  sinusoid  one  was  varied  from  0.1  to  100  and  the 
power  of  sinusoid  two  was  held  constant  at  0.1.  The  threshold  power  was 
1.  The  major  weaknesses  of  the  APS2  are  apparent  from  this  figure: 

(A)  Because  the  gain  of  a  DFT  bandpass  filter  is  not  unity  for  off-bin 
signals,  the  gain  at  the  frequency  of  sinusoid  one  remained  unity  when 
the  power  of  the  sinusoid  v/as  between  1  and  2,  even  though  the  threshold 
power  was  1.  (B)  Because  of  poor  filter  design,  when  the  gain  at 

sinusoid  one's  frequency  did  change,  it  first  increased  and  then  slowly 
dropped  towards  zero,  rather  than  immediately  becoming  zero.  (The  only 
reason  that  the  gain  drops  as  the  signal  power  increases  is  that  leakage 
causes  more  and  more  of  the  poles  to  be  removed.)  (C)  The  gain  at  the 
frequency  of  the  second  sinusoid  dropped  to  zero  even  though  the  power 
at  that  frequency  was  always  less  than  the  threshold  power.  This  is 
the  result  of  leakage  from  sinusoid  one. 

The  effects  of  leakage  can  also  be  seen  in  Figure  4-10  which  shows 
the. transfer  function  of  the  APS2  for  several  powers  of  signal  one. 
Figure  4-10  show  that  the  APS2  is  an  all  pass  filter  when  the  power  of 
signal  one  is  0.1.  When  the  power  of  signal  one  is  2.5,  the  APS2 
develops  a  notch  near  the  frequency  of  signal  one.  However,  sidelobe 
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Transfer  function  of  APS2  for  various  input  powers. 


effects  cause  the  gain  at  the  frequency  of  signal  one  to  be  greater  than 
one.  Vlhen  the  power  of  signal  one  is  100  the  APS2  notches  a  wide  range 
of  frequencies,  including  the  frequency  of  signal  two.  This  wide  notch 
is  caused  by  leakage  from  signal  one. 
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5  --  APS3:  AN  LALE-DFT  BASED  ADAPTIVE  POWER  SEPARATOR 

Introduction 

This  chapter  proposes  and  examines  an  APS  which  uses  the  LALE  and 
a  DFT  for  spectral  analysis,  and  a  DPT  based  filter  design  technique  . 
called  frequency  sampling  [Rab]  to  design  bandstop  filters.  This 
adaptive  power  separator  will  be  called  the  APS3  to  distinguish  it  from 
the  purely  DFT  based  APS2  that  was  analyzed  in  Chapter  4.  The  reason 
for  using  a  LALE  and  a  DFT  for  spectral  estimation  is  that  the  LALE 
reduces  leakage  effects  and  is  not  sensitive  to  the  frequency  of  the 
signals  (i.e.  the  signals  do  not  have  to  be  'on-bin').  The  advantage 
of  the  frequency  sampling  filter  design  technique  is  that  it  is  simple, 
fast,  and  by  increasing  the  filter  length  can  create  arbitrarily  good 
filters.  Thus  the  APS3  improves  on  the  two  weaknesses  of  the  APS 2: 
poor  power  spectral  estimation  and  poor  filter  design. 

The  APS3  is  based  on  the  same  simple  idea  as  the  APS2.  Estimate 
the  power  spectrum  of  the  input  and  then  design  a  filter  that  passes 
the  frequencies  which  have  low  power  signals  and  rejects  the  frequencies 
which  have  high  power  signals.  Figure  5-1  shows  a  block  diagram  of  the 
APS3.  The  theory  of  operation  is:  1)  the  LALE  forms  a  filter  with  a 
gain  that  is  related  to  the  signal -to-noise  ratio  (SNR).  The  larger 
'  the  signal  power  the  larger  the  SNR  and  consequently  the  larger  the 
filter  gain.  2)  The  DFT  of  the  weight  vector  of  the  LALE  measures  the 
gain  of  the  LALE,  which  is  a  modified  power  spectral  estimate.  3)  The 
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Block  Diagram  of  the  APS3. 
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modified  power  spectral  estimate  is  used  to  create  a  desired  filter 
characteristic,  which  is  then  realized  as  a  practical  filter.  4)  The 
resulting  filter  is  used  to  filter  the  input,  and  by  design  the  output 
of  this  filter  will  contain  only  the  weak  signals. 

The  next  few  sections  are  devoted  to  analyzing  the  modified  power 
spectral  estimate  that  is  derived  from  the  LALE.  We  show  that  this 
modified  power  spectral  estimate  has  less  leakage  than  a  DFT  based  power 
spectral  estimate.  Thus  the  LALE  power  spectral  estimate  can  detect 
small  signals  that  are  near  powerful  off-bin  signals  better  than  a 
DFT  power  spectral  estimate.  (Mote:  the  use  of  windows  with  the  DFT 
will  reduce  the  effects  of  leakage,  but  very  powerful  off-bin  signals 
will  still  swamp  small  signals.)  Thus  the  LALE  power  spectral  estimate 
can  be  used  to  obtain  arbitrarily  good  switching  characteristics,  so 
the  APS3  can  have  a  lower  power  threshold  p-j  that  is  arbitrarily  close 
to  the  upper  power  threshold  p2  (see  the  descriptions  of  tolerance 
scheme  one  and  two  which  are  given  by  Eq.  3-1  and  3-2). 

The  last  sections  of  this  chapter  analyze  the  filter  design 
technique  used  in  the  APS3.  Vie  prove  that  the  filter  design  technique 
used  in  the  APS3  can  design  arbitrarily  good  (in  an  integral  squared 
error  sense)  filters.  Thus  we  show  that  the  APS3  fits  gain  tolerance 
scheme  one  (Eq.  3-1 ) . 

LAL E-DFT  Based  Power  Spectral  Estimate 

«.  In  this  section  we  will  use  the  analysis  of  the  mean  LALE  weight 

vector,  which  was  presented  in  Chapter  2,  to  show  that  the  DFT  of  the 
LALE  weight  vector  is  a  modified  power  spectral  estimate.  The  perform¬ 
ance  of  the  LALE-DFT  power  spectral  estimate  will  be  analyzed  and 
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compared  to  the  conventional  Welch  periodogram  power  spectral  estimate. 
Again,  the  analysis  in  this  section  follows  Treichler's  [Tre]  analysis 
with  only  small  mod  if ications  to  account  for  the  use  of  a  LALE  as 
opposed  to  an  ALE. 

The  discrete  power  spectrum  of  a  signal  x(k)  may  be  found  by  taking 
the  OFT  of  n  samples  or  the  autocorrelation  function  of  x(k).  If  the 
autocorrelation  of  x(k)  is  r(t),  then  samples  of  the  autocorrelation 
function  can  be  grouped  to  form  the  vector:  4 

T  =  C r( A ) *  r(A+l),  ....  r(A+n-l )]T 

The  discrete  power  spectrum,  S,  corresponding  to  this  autocorrelation 
vector  is: 

S  =  F-T  (5-1) 


where  F  is  an  m  by  n  matrix  formed  from  the  basis  vectors  of  the  DFT . 

The  elements  of  the  F  matrix  are: 

[F]i  £  =  Exp{j2Tr(~  +  <j>) }  (5-2) 

where  the  starting  phase,  <p,  is  arbitrarily  set  to  correspond  to  the 
decorrelation  delay  A- 

The  ALE  spectral  estimate,  which  is  derived  from  the  DFT  of  the 
expected  LALE  weight  vector,  can  be  written  as: 


'LALE 


=  F  ■  E { W) 
n-1 

-  f  Z  T - 

1=0  °  +  Y  +  Xi 


[CL 


(5-3) 
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where  r  is  the  autocorrelation  vector  of  a  sinusoid 


C 


p  is  the  power  of  the  sinusoid 
Then  for  a  sinusoid  at  frequency  u». 


Itu],! 


SinfJ^  -  uiT ) }  Sin{£<2^  +  WT)} 

- - - - - -  +  — - T - t - - - -  - 

nSin{^(2^  -  WT)}  nSin{|-(27;i  +  tuT)> 


(5-7) 


For  comparison,  the  Welch  spectral  estimate  is: 

K 

S  =  £  5  P.-IL  .  (5-9) 

i=l  ^  11 

Thus  the  LALE-DFT  produces  an  amplitude  distorted  power  spectral 
estimate  when  compared  to  the  Welch  spectral  estimate.  However,  the 
'  amplitude  distortion  can  be  compensated  for  in  the  APS3  by  adjusting 
threshold  power  p^.  Since  the  amplitude  of  the  LALE-DFT  power 
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spectra]  estimate  is  a  nonotonic  function  of  the  signal  power,  an 


such  that  if  the  amplitude  of  the  LALE-DFT  power  spectral  estimate  is 
greater  than  a^.  then  the  power  of  the  input  signal  must  be  greater  than 
Pt-  Thus  the  LALE-DFT  power  spectral  estimate  can  be  used  by  the  APS3 
to  distinguish  between  weak  signals  and  strong  signals. 

Compared  to  the  Welch  spectrum,  the  LALE-DFT  spectrum  provides 
better  detectability  of  relatively  weak  signals  that  are  close  in 
frequency  to  powerful  signals.  This  advantageous  behavior  is  a  result 
of  the  biased  amplitude  in  the  LALE-DFT  spectrum.  Recall  that  the 
LALE-DFT  spectrum  is  composed  of  a  sum  of  scaled  shape  functions 
(Eq.  5-8).  These  shape  functions,  U,  have  a  maximum  amplitude  of  1  and 
sidelobes  whose  amplitudes  are  inversely  proportional  to  the  difference 
between  the  frequency  of  the  sidelobe  and  the  mainlobe.  The  maximum 
scale  factor  is  unity,  so  the  largest  possible  peak  in  the  LALE-DFT 
spectrum  has  an  amplitude  of  unity.  Therefore,  since  the  largest  side- 
lobe  of  the  shape  function  is  typically  less  than  0.25  for  a  reasonably 
sized  DFT,  the  largest  sidelobe  in  the  LALE-DFT  spectrum  is  smaller 
than  0.25.  If  y  is  appropriately  chosen  so  that  a^,  the  amplitude 
threshold,  is  greater  than  0.25,  then  no  sidelobe  can  exceed  at  thus  no 
sidelobe  can  be  falsely  treated  as  a  powerful  signal.  (However,  side- 
lobes  can  add  and  thereby  cause  a  weak  signal  to  appear  strong  or  a 
strong  signal  to  appear  weak.  This  effect  is  relatively  minor  and  only 
of  concern  when  there  are  closely  spaced  input  signals,  consequently 
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the  superposition  of  sidelobes  can  usually  be  ignored  if  all  of  the 
input  signals  are  separated  by  at  least  1.5/nT  Hz.)  In  contrast  the 
Welch  spectrum  is  also  composed  of  scaled  shape  functions,  but  the 
scaling  factor  is  proportional  to  the  signal  power.  Thus  a  powerful 
signal  will  have  large  sidelobes  which  would  be  falsely  treated  as 
powerful  signals. 

Figure  5-2  illustrates  the  difference  between  the  LALE-DFT  spectrum 
and  the  Welch  spectrum.  The  input  consisted  of  the  sum  of  two  sinusoids. 
One  sinusoid  was  at  25%  of  the  Nyquist  frequency  with  a  power  of  1 ,  and 
the  other  was  at  32.8125%  of  the  Nyquist  frequency  with  a  power  of  10. 
These  spectral  estimates  were  computed  with  a  256  point  DFT  of  a  32 
weight  LALE  (224  zeros  are  appended  to  the  32  values  of  the  LALE  weights 
to  generate  a  256  point  vector  for  the  DFT)  and  a  256  point  DFT  of  a 
64  point  autocorrelation  function,  since  a  64  point  autocorrelation 
function  is  as  easy  to  compute  as  a  32  weight  LALE.  (Note:  the  DFT  of 
a  64  point  autocorrelation  function  has  better  frequency  resolution  than 
the  DFT  of  a  32  point  autocorrelation  function,  thus  this  comparison  is 
not  biased  in  favor  of  the  LALE.)  The  LALE  power  spectrum  clearly  shows 
the  peak  due  to  the  smaller  signal,  whereas  the  smaller  signal  is 
swamped  in  the  Welch  spectrum.  Thus  for  any  fixed  computational  com¬ 
plexity,  the  LALE-DFT  power  spectral  estimate  has  less  leakage  than  the 
Welch  spectrum. 

We  have  established  that  the  LALE-DFT  power  spectrum  can  be  used 
».  to  determine  the  frequencies  at  which  the  signal  power  exceeds  a  pre¬ 
determined  threshold.  However,  if  a  simple  threshold  device  is  used  to 
distinguish  between  low  power  and  high  power  signals,  noise  in  the 
spectral  estimate  could  cause  serious  problems.  For  example,  if  the 
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Welch  spectrum 


LALE  spectrum 

Power  of  signal  at  frequency  fj  =  1.0 
Power  of  signal  at  frequency  f^  =  10.0 


5-2.  Comparison  of  Welch  and  LALE  power  spectra. 


( 


expected  spectrum  is  near  the  threshold  value  then  noise  could  cause 
the  threshold  device  to  switch  repeatedly  as  the  actual  spectral 
estimate  varied  about  the  expected  value.  Switching  by  the  threshold 
devices  would  cause  the  APS3  to  redesign  its  bandstop  filter,  and  this 
would,  in  turn,  modulate  the  output  which  is  often  undesirable.  Thus 
v/e  need  to  analyze  the  noise  in  the  power  spectral  estimate. 

Variance  of  the  LALE  weight  vector 

The  variance  of  the  LALE-DFT  power  spectral  estimate  depends  on  the 
variance  of  the  LALE  weight  vector,  so  in  this  section  we  analyze  the 
variance  of  the  LALE  weight  vector.  As  in  Chapter  2,  this  analysis 
follows  Treichler's  [Tre]  analysis  with  only  those  modifications  required 
to  account  for  the  change  in  adaptive  algorithm. 

Define  the  weight  vector  noise  to  be; 

( 

AVl(k)  =  W(k)  -  E{W( k) }  (5-11) 

then  the  weight  vector  covariance  is: 

E{AW( k) AWT( k) }  =  n(k)  .  (5-12) 

The  evolution  of  AW(k)  is  complicated,  so  we  study  instead  the  evolution 
of  V(k)  where: 

v(k+l )  =  [vI-2yR(k)]V(k)  +  2y[d(k)X(k)-Ps]  .  (5-13) 

The  evolution  of  V{k)  closely  models  the  average  behavior  of  AW(k), 
especially  if  y  is  small  and  the  LALE  is  near  convergence.  Since  V(k) 
closely  approximates  AW(k),  the  converged  covariance  of  V  will  approxi¬ 
mates  .  W(k),  the  converged  covariance  of  V  will  approximate  the  converg¬ 
ed  covariance  of  AW. 

The  covariance  of  V  can  be  written  as: 
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Cov { V ( k+ 1 ) }  =  E{V(k+l )VT( k+1 } } 

=  [vI-2yR( k)]Ccv { V ( k) }  [vI-2yR(k}] 

k 

+  4u2[E{B(k,0)}  +  5>([vI-2pR( k-i  )]n 
i=l 

+  E{BT(k,i)}[vI-2pR(k,i)]1')]  .  (5-14) 

where  the  matrix  B{k,i)  is  defined  as: 

B(k, i)  =  [d(k)X(k)-Ps]  [d(k-i)X(k-i)-Ps]T  .  .(5-15) 

Equation  5-15  is  a  function  of  the  specific  input  signals  and  very 

2 

complicated.  However,  if  the  input  consists  of  white  noise  of  power  o 
and  a  sinusoid  of  power  p  and  angular  frequency  u  then  Treichler  [Tre] 
showed  that: 

2 

[£{B(k,i)}]^m  =  f  Cos[ulUHa-21)3 

+  pCos[wT(£-m)]a2o(£-m+i) 

2 

+  pCo  s [wT ( £-m+i ) ] o  6 ( i ) 

+  a46(i)5U-m+i) 

By  transforming  Eq.  5-14  with  the  modal  matrix  Q>  and  defining 
V(k)  =  QV'(k),  n  =  Qn'QT,  and  assuming  that  u  and  Y  are  small  enough 
so  that  [vI-2pR]n+A  r  I,  then: 

E {V ' ( k+1 ) V 1 T( k+1 ) }  =  [vI-2u(o2+A)]E{V (k)V,T(k) } 

+  [vI-2p(o2+A}] 

+  4p2(2a2A+o4I)  .  (5-16) 
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Both  the  transition  matrix  and  the  driving  term  of  Eq.  5-16  are 
diagonal,  so  the  converged  value  of  n‘  is  diagonal.  Thus  Eq.  5-16  can 
be  solved  for  the  converged  value  of  ft': 


ya2(2Xi+o2) 

(Y+a  +Xj)  -  y(Y+a  +X^)‘ 


2  xi  ~  y 

Z  yc>  (1  +• - - ) 


Y  +  0  +  Y ; 


(5-17) 


Equation  5-17  can  be  solved  for  n(°°)  by  reversing  the  coordinate 
transformation: 


n(«)  =  va2(i  +  £ - h — -QiQl)  * 

i-0  y  +  a  + 


(5-18) 


Equation  5-18  is  valid  if  the  input  consists  of  multiple  sinusoids  in 

white  noise.  As  a  check,  if  the  input  has  low  SNR,  i.e.  o'  » 

2 

then  n(°°)  ~  uo  I  which  is  the  result  derived  in  [Wid]  for  stochastic 
inputs. 

Variance  of  the  LALE-DFT  power  spectral  estimate 

Denote  the  actual  LALE-DFT  spectral  estimate  as  and  the 

a 

mean  spectral  estimate  as  _  F-E{W}.  Using  these  definitions  the 

covariance  of  is: 


C0V{SLALE}  =  E^SLALE'SLALE^SLALE'SLALE^  } 


=  F  E{(W-E{W})(W-E(W})T}FT 


(5-19) 
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From  Eq.  5-18, 


where: 


n  =  £{(w-e{w})(w-ew;t} 


ua  (I+Re) 


(5-20). 


< 

V 


n_1  A- 

£ - jJ— 

i=0  y  +  o  +  A. 


q.qT 


=  coherent  input  equivalent 
autocorrelation  matrix 


Substituting  Eq.  5-20  into  5-19  yields: 


C0V{\ALE}  =  Ma2(nI+FReFT) 


(5-21) 


jr 

\ 


The  variance  of  the  spectral  estimate  is  the  main  diagonal  of 
C°v t^iALE^ '  Solving  for  the  variance  yields: 


[Var{SLALE}]m 


n-1  n-1 

=  Z  Z  [nL  .  Exp(2Trjm— -^) 
i=0  1=0  1,4 


=  pa 


•  n-1 

1  +  Z  nzir1  r  e(  t)  Exp(  2tt  jm^-) 

t~- ( n-1 )  n  e  n 


where 


(5-22) 


r  (t)  =  the  coherent  component  equivalent 

e  autocorrelation  function 
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Equation  5-22  shows  that  the  variance  of  the  spectral  estimate  is 
a  constant  plus  a  curve  which  is  determined  by  the  DFT  of  a  Bartlett- 
windowed  version  of  the  equivalent  coherent  autocorrelation  function 
re(t)-  An  important  property  of  the  variance  is  that  it  is  proportional 
to  y.  Thus  by  decreasing  y,  the  variance  of  the  spectral  estimate  can 
be  reduced  to  any  desired  value. 

A  simple  example  serves  to  illustrate  the  properties  of  the 

4  mn 

variance.  Suppose  that  the  input  is  a  sinusoid  of  frequency  — - 

nr 

(an  on-bin  signal)  and  power  p,  then: 


rvar'W]m 


nycr  1  + 


n-1 

_ P _ _  sp  n  - 

y  +  a2  +  £p  t=-(n-l )  n 


Cos(2u-y-)  Exp(27Tjm£) 


nyo  nixing  or  n-mQ 


2  5  ESNR 

nyo  (1  +  - = - )  m=mn  or  n-mn 

1  +  ~  ESNR  0  0 


(5-23A) 


(5-23B) 


For  the  design  of  the  APS3,  Eqs.  5-23A  and  5-23B  can  be  reduced  to  one 
simple  design  rule: 


2  2 
nyo  <  Var(sLALE>  £  2nya 


(5-24) 


In  summary,  we  have  found  that  if  the  input  could  be  modeled  as  a 
sum  of  sinusoids  in  white  noise  then  the  LALE-DFT  spectrum  has  the  same 


shape  as  the  ideal  power  spectrum,  but  that  the  power  estimate  is 
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biased  (Eq.  5-15).  Furthermore,  we  have  shown  that  the  variance  of  the 
spectral  estimate  can  be  bounded  by  a  simple  expression  (Eq.  5-24).  In 
the  next  section  we  will  use  these  results  to  show  that  the  LALE-DFT 
power  spectrum,  in  conjunction  with  a  hysteresis  threshold  device,  fits 
the  power  switching  requiranents  of  both  gain  tolerance  schemes  that 
were  presented  in  Chapter  2. 


How  the  LALE-DFT  fits  the  gain  tolerance  schemes 

In  this  section  we  will  show  how  the  LALE-DFT  power  spectral 

estimate  can  be  used  in  the  APS3,  and  how  the  APS3  will  meet  the  switch¬ 

ing  requiranents  of  both  gain  tolerance  schemes.  This  is  an  importnat 
step  towards  a  useful  adaptive  power  separator. 

Both  gain  tolerance  schemes  have  two  power  thresholds,  p.j  and  p2-  - 
A  switching  device  that  obeys  the  tolerance  schemes  should  indicate 
"pass"  if  the  input  signal  power  is  less  than  Pp  "stop"  if  the  signal 
power  is  greater  than  p2  and  either  "pass"  or  '’stop"  otherwise.  From 
Eq.  5-10  we  find  that  the  amplitude  of  the  LALE-DFT  spectrum  is  a 
monotonic  function  of  the  input  signal  power.  Thus  a  unique  spectrum 
amplitude  is  associated  with  each  input  power.  Assume  that  an  amplitude 

threshold  at  is  found  (using  Eq.  5-18)  which  corresponds  to  an  input 

(Pi+P2) 

power  of  — ^ .  Then  if  the  amplitude  of  the  LALE-DFT  spectrum  is 

(Pl+P2) 

greater  than  a^  the  input  power  must  be  greater  than - ^ - .  Thus  if 

n  threshold  devices  with  equations: 
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are  used  on  the  output  of  the  LALE-DFT  they  will  indicate  which  ( 

frequencies  should  be  passed  and  which  should  be  stopped. 

Noise  in  the  spectral  estimate  will  adversely  affect  the  threshold 
devices.  For  example,  because  of  noise  the  amplitude  of  the  spectral 
estimate  may  be  less  than  a^  even  though  the  signal  power  is  greater 
(pl+P2) 

than  — - .  This  will  cause  the  threshold  device  to  output  a  "pass.'’ 

when  it  should  indicate  "stop."  Thus  with  a  noisy  spectral  estimate 
there  is  a  possibility  of  a  false  output  from  the  threshold  devices. 

Since  the  statistics  of  the  spectral  noise  are  known,  these  effects  can 
be  quantified. 

In  accordance  with  the  gain  tolerance  schemes,  v/e  say  there  is  an 
error  (false  output)  when  the  threshold  devices  indicate  "pass"  even 
though  the  input  power  is  greater  than  p^  or  when  the  thresholc  devices  ^ 

indicate  "stop"  even  though  the  input  power  is  less  than  p^ .  The  proba¬ 
bility  of  an  error  is: 


( 
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Thus,  provided  that  f  ,  the  probability  of  an  error  is  finite 
and  can  be  made  as  small  as  desired  by  making  p  small.  Therefore  the 
LALE-DFT  spectral  estimator  fits  both  tolerance  schemes  in  a  probability 
of  an  error  have  to  be  balanced  against  convergence  time. 

Earlier  we  discussed  the  possibility  that  a  signal  with  power  near 
the  threshold  power  could  cause  the  output  of  the  threshold  devices  to 
repeatedly  change  state.  This  could  be  undesirable  since  it  would 
cause  the  APS3  to  repeatedly  redesign  its  bandstop  filter  which  would, 
in  turn,  modulate  the  signals  that  pass  through  the  filter.  A  measure 
of  this  behavior  is  the  expected  design  cycle  time,  where  the  design 
cycle  time  is  the  time  it  takes  for  the  output  of  a  threshold  device  to 
change  state  and  then  return  to  the  original  state.  If  the  design  cycle 
time  is  large  then  the  filter  is  infrequently  redesigned  and  will  not 
continually  modulate  the  signals.  Thus  we  will  analyze  the  APS3 
threshold  devices  to  find  the  design  cycle  time  T^. 

The  output  of  a  threshold  device  can  be  viewed  as  a  Markov  process. 
The  output  has  two  states,  "pass"  and  "stop,”  and  for  each  input  power 
there  are  probabilities  associated  with  each  state.  We  define  Pp  as 
the  probability  that  the  output  is  a  "pass"  and  P$  as  the  probability 
that  the  output  is  a  "stop."  The  state  diagram  for  the  threshold 
device  is  shown  in  Figure  5-3.  The  expected  design  cycle  time  can  be 
computed  [ Ci n] : 


E {Tc >  =  E{time  from  pass  to  stop}  +  E{time  from  stop  to  pass} 


-L  +  -L 

PP  Ps 


(5-29) 


c 
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If  the  signal  power  is  exactly  the  threshold  power,  Pp  =  P$  ■  0.5,  and 

so  E{Tc}  =  4  iterations.  If  this  time  is  too  short,  a  more  complicated 

threshold  device  will  have  to  be  used  because  there  is  no  way  to  control 

Pp  and  P$  with  a  simple  threshold  device. 

A  threshold  device  with  hysteresis  will  allow  the  designer  to 

control  Pp  and  P$.  By  making  Pg  or  Pp  small,  E{Tc >  can  be  made  as 

large  as  desired.  To  see  how  a  hysteresis  device  provides  control  over 

P,  and  P  we  define: 
s  p 


<|>(k)  =  output  of  the  threshold  device  at  time  k 


i 


SLALE 


input  to  the  threshold  device 


a^  =  lower  threshold 


ah  =  higher  threshold 


Using  these  definitions  the  law  for  the  hysteresis  threshold  device  is: 


$(k+l ) 


pass" 

i  <  a^ 

pass" 

al  <  1  <  ah 

and 

4>(k)  =  "pass" 

stop" 

ah^ 

stop" 

a£  <  i  <  flh 

and 

4>(k)  =  "stop" 

(5-30) 


The  probabilities  of  changing  state  are  conditioned  on  the  current 
S,tate,  $(k): 
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ps|p  s  p^tk+1)  =  ”stop"|o(k)  =  "pass”}  =  P{i>ah} 

P_l„  =  p($(k+l)  =  “pass"  |${k)  =  "pass")  =  P{i<a.  } 

P(P  n 

Ps|s  =  p{*<k+1>  =  “stop"  |<j>(k)  =  "stop"}  =  P{i>ai} 

» 

Pp|s  *  P{<t>(k+1 )  =  "pass" |<t»(k)  =  "stop"}  *  Ptffaj}  .  .  (5-31) 

The  state  diagram  for  the  hysteresis  device  is  shown  in  Figure  5-4. 

For  this  device  the  expected  cycle  time  is: 

E(TJ  »  ,-!-  *  jJ—  .  (5-32) 

C  Ps|p  Pp|* 


THEOREM  5-1:  If  a^  f  a^  then  E {Tc >  can  be  made  as  large  as  desired. 

PROOF:  If  either  P  ,  or  P  .  can  be  made  as  small  as  desired  then 

s  |p  pjs 

E(TC)  can  be  made  as  large  as  desired,  so  we  will  ;how  that  either 

P  i  or  P  .  can  be  made  as  small  as  desired, 
sjp  p|s 

Let  6  «  Max(E{i}-a^,ah-E{i}).  Assume  for  convenience  that 
5  =  E{i}-a^,  then: 


P(i£a^} 

P{E{i}-i>6} 


Equation  5-24  shows  that  Var{S^Af can  be  made  as  small  as 
desired  by  making  the  value  of  y  small.  Therefore  Ps|p  can 
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made  as  small  as  desired,  and  consequently  E{Tc>  can  be  as  large 
as  desired  since: 

E(T)  >  5^—  ♦  1  . 

C  Ps|p 

An  analagous  demonstrations  can  be  used  to  show  that  Pp|s  can  be 
made  as  small  as  desired  by  defining  d  *  a^  -  E{i). 

This  completes  the  proof  of  theorem  5-1. 

The  expected  design  cycle  time  can  be  bounded  using  the  same  type 
of  arguments  that  were  presented  in  the  proof  of  theorem  5-1  .  A  lower 
bound  on  E {T£ }  is: 


E{T  }  > 


c  - 


ah  -  a£ 

2VAR{W 


ah  '  *1 
2 

4npo 


(5-33) 


When  a  hysteresis  device  is  used  in  the  APS3,  the  probability  of 
an  incorrect  state  or  error  is: 

P{error)  <  P((i>ah  and  P<P-| )  or  (i<a£  and  p>p2)> 

<  P{|t-E(U|>mn(.h-ainjn. 


Var(W 


Min(ah-a  aM  *a.) 
n  min  max  L  • 


(5-34) 
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where 


c 


a 


min 


2?1 
+  Y  + 


Frequency  sampling  filter  design 


In  the  previous  sections  we  have  analyzed  the  power  spectral 
estimation  technique  used  by  the  APS3.  We  have  shown  that  the  DFT  of 
the  LALE  weight  vector  in  conjunction  with  hysteresis  threshold  devices 
can  indicate  with  any  desired  accuracy  which  frequencies  should  be 
blocked  and  which  should  be  passed.  In  the  following  section  we  will 
show  that  the  frequency  sampling  filter  design  technique  can  use  the 
output  of  the  threshold  devices  to  design  a  filter  that  will  block,  to 
any  desired  accuracy,  the  high  power  signals  and  pass,  with  nearly  unit 
gain,  the  low  power  signals. 

One  technique  for  designing  finite  impulse  response  (FIR)  filters 
is  the  frequency  sampling  technique  [Opp,  Rab,  Rab2],  Frequency  sampl¬ 
ing  uses  an  inverse  OFT  to  obtain  a  weight  vector  from  samples  of  the 
ideal  frequency  response.  This  technique  is  fast,  computationally 


r 

V 
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efficient  and  not  iterative.  Since  the  APS 3  will  be  a  real-time  signal 
processor  speed  is  important.  Also,  frequency  sampling  can  be  used  to 
design  filters  with  linear  phase  which  may  be  desirable  in  communication 
systems.  Therefore  the  frequency  sampling  technique  is  used  in  the  APS3 
to  design  the  bandstop  filter  that  rejects  the  strong  input  signals. 

In  the  APS3  the  ideal  frequency  response  is  derived  from  the  DFT 
of  the  LALE  weight  vector  n  threshold  devices  process  the  outputs  of  an 
n  point  DFT  to  produce  samples  of  the  ideal  frequency  response.  The 
frequency  sampling  technique  uses  the  following  equation  (inverse  DFT) 
to  obtain  a  weight  vector  from  the  sampled  frequency  response: 

n-1 

[ML  =  ~  £  [H]»Exp(2irj  ^-)  i  *  0,  1,  ....  n-1  (5-35) 

1  n  1=0  *■  n 

where 

W  =  the  weight  vector  of  the  FIR  filter 
=  the  pulse  response  of  the  FIR  filter 
H  =  the  vector  of  gain  specifications 

*  samples  of  the  ideal  frequency  response 
[H]^  *  ideal  gain  at  frequency 

n  »  number  of  weights  in  FIR  filter 
=  number  of  frequency  samples 
T  =  sampling  period  of  FIR  filter 

Since  we  stipulate  that  W  must  be  real,  the  following  symmetry 
condition  has  to  be  imposed  on  H: 

[H]£  =  Conjugate([H]n_^)  .  (5-36) 


The  transfer  function  of  an  FIR  filter  designed  by  Eq.  5-35  is  [Opp3: 


H(Z)  = 


1  -  Z' 


t‘ 0  1  -  ExpUZti^Z"1 


(5-37) 


The  frequency  response  of  the  FIR  filter  is  found  by  substituting 
Z  =*  Exp(jti))  into  Eq.  5-37: 

.  ,  .  n-1  »  ,  Sin[y(u>-2ir|)3 

H(ejw)  =  I  Exp{-ja^-)  T,  (-1  r[H)£Exp(  ji4)  - f- - 5— 

n  2  £o  t  n  sin[i(u»-2TrJ)3 


(5-38) 


substituting  u>  ■  tt—  into  Eq.  5-38  shows  that: 


jiA 

[H(e  ")| 


*  |[H]£|  l  =  0,  1 . n-1 


Thus  the  gain  of  the  FIR  filter  has  the  correct  magnitude  at  the  sample 
frequencies  (integer  multiples  of  Af).  What  about  the  phase?  Eq.  5-37 
shows  that  the  frequency  response  consists  of  a  linear  phase  term, 

Exp(-jury^),  and  a  complex  term.  One  of  the  design  requirements  is 
that  the  FIR  filter  have  linear  phase.  This  will  be  true  if  the  complex 
term  is  purely  real  which  requires  that: 


[G]£  e  n 


(5-39) 


s  where  [G]£  a  rea1  value,  the  magnitude  of  the  ideal  frequency 
response.  Also,  since  H  has  to  obey  Eq.  5-36: 

[G]»  =  -[G]n  , 


(5-40) 


98 


If  H  has  the  form  specified  in  Eqs.  5-39  and  5-40,  then  the  FIR  filter 
will  have  a  real  weight  vector,  linear  phase  and  the  same  magnitude  as 
the  ideal  filter  at  the  specification  frequencies  (integer  multiples  of 

4j-Hz).  However,  there  will  be  zeros  in  any  passband,  as  can  be  seen 
by  rewriting  Eq.  5-38  as: 


H(eJ“) 


n-1 


L  (-I)^CGL 
Z=0  4 


Sinfflu-Znft] 

nSin[J-(w-2ir^)] 


(5-41) 


This  shows  that  the  magnitude  of  H(Z)  is  the  sum  of  n  bandpass  filters. 

t 

The  peak  amplitude  of  the  filters  is  [G]^  at  frequency  irj-.  In  a  pass- 
band,  where  [G]^  ^  =  [G]^  =  1,  the  sign  alternations  in  H(eju)  caused 

by  the  term  (-1)  will  cause  the  gain  to  be  zero  somewhere  between 
frequency  and  This  follows  from  the  intermediate  value  theorem  ( 

[Ros].  Since  this  behavior  is  undesirable,  an  additional  condition  is 
imposed: 

CG]£(-1)^>0  .  (5-42) 

With  the  added  restriction  of  Eq.  5-42  the  magnitude  of  H(Z)  will 
approximate  quite  closely  the  magnitude  of  the  ideal  frequency  response. 

Figure  5-5  shows  the  magnitude  of  H(Z)  for  a  typical  bands top  filter 
designed  by  frequency  sampling.  The  frequency  response  has  ripples 
which  increase  in  magnitude  near  discontinuities  in  the  derivative  of 
'  the  ideal  frequency  response. 

The  magnitude  of  the  ripples  in  the  actual  frequency  response  can 
be  reduced  by  smoothing  the  ideal  frequency  response.  A  common  technique 

C 
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for  smoothing  the  ideal  response  is  to  vary  a  few  of  the  ideal  gains 
near  a  discontinuity.  Usually  the  gains  are  varied  until  the  ripple  is 
minimized  [Opp,  Rab,  Rab2].  This  procedure  is  complicated  and  nay  be 
too  slow  for  real-time  filter  design.  Thus  the  APS3  uses  fixed  smooth¬ 
ing  of  the  ideal  filter  characteristic.  Wherever  there  is  a  transition 
from  a  desired  gain  of  one  to  a  desired  gain  of  zero,  the  APS3  will 
alter  the  gain  of  one  to  a  gain  of  0.4  (0.4  was  chosen  by  examining  a' 
table  of  optimal  values  [Rab2]  and  selecting  a  compromise  value). 

Some  reduction  in  the  complexity  of  the  APS3  can  be  obtained  by 
using  the  symmetry  of  Eqs.  5-39  and  5-40  and  the  fact  that: 


if 


W  »  inverse  DFT(H) 


(5-43A 


then 


U'  =  inverse  DFT(H' ) 


(5-433 


where 


[w']i  •  [W] 


i  +  j  mod  n 


on.  -  (-D1’[H]i 

Figure  5-6  shows  a  block  diagram  for  one  possible  implementation  of  the 
filter  design  portion  of  the  APS 3.  Notice  how  symmetry  eliminates  the 
need  for  half  of  the  threshold  devices  and  half  of  the  smoothers.  Alsc 
the  symmetry  conditions  required  for  linear  phase  dictate  that  the  n/2 
input. must  be  zero  since  [G]  *  -[G]  from  Eq.  5-40.  Finally,  the 


output  of  the  inverse  OFT  is  'rotated',  which  by  Eqs .  5-43A  and  B  has 
the  same  effect  as  Eq.  5-42. 

This  completes  the  derivation  of  the  frequency  sampling  technique 
for  filter  design.  In  the  next  section  we  will  sho w  that  the  frequency 
response  of  filters  designed  by  the  APS3  can  be  made  as  nearly  ideal 
as  desired. 


Quality  of  frequency  sampling  designed  filters 

The  quality  of  a  filter  can  be  defined  in  many  ways-  One  measure 
of  the  quality  of  a  filter  is  defined  as  the  integral  of  the  square  of 
the  difference  between  the  ideal  transfer  function  and  the  actual 
transfer  function.  The  following  theorem  uses  the  integral  squared 
error  [Trt,  Oga]  criterium  to  characterize  a  frequency  sampling  filter. 

THEOREM  5-2:  The  integral  squared  error  between  the  transfer  function 
of  a  filter  designed  by  frequency  sampling  and  the  ideal  transfer 
function  can  be  reduced  to  any  desired  value  by  making  the  filter 
length  large  enough. 

A  useful  corollary  to  Theorem  5-2  is: 

COROLLARY  5-2:  Given  a  signal  with  a  finite  power  spectral  density 
and  a  frequency  sampling  designed  FIR  filter,  the  output  power  due  to 
the  difference  between  the  actual  transfer  function  and  the  ideal 
transfer  function  can  be  made  as  small  as  desired. 

* 

Before  we  prove  Theorem  5-2  we  will  prove  two  lemmas. 


LEMMA  1:  The  impulse  response  of  a  frequency  sampling  designed  FIR 
filter  is  an  aliased  or  folded  version  of  the  impulse  response  of  the 
ideal  filter  [Rab2]. 

PROOF:  (The  technique  used  in  this  proof  is  due  to  Papoulis  [Pap].) 

The  transfer  function  H ( u>) ,  of  a  digital  filter  is  periodic  with  a 
period  of  y.  This  implies  that  H(u>)  can  be  expanded  in  a  Fourier 


senes: 


H(cj)  =  E  c^Exp(-j2mkTu) 
k=-“ 


where  the  c^,  are  the  Fourier  coefficients  and  also  the  pulse  response 
of  the  ideal  filter. 

Frequency  sampling  uses  samples  of  the  ideal  transfer  function. 
Let  n  be  an  arbitrary  integer  and  define: 


°1  "  nT 


then  the  samples  of  the  ideal  transfer  function  used  by  frequency¬ 


sampling  are: 


HCmu-j )  m  =  0,  1 ,  . .. ,  n-1 


In  terms  of  the  Fourier  coefficients  these  samples  are: 


H(mo.)  =  E  Cj^Expt-jETrkTfTK^) 

k=-°° 


E  c.Exp(-j2-kJ)  m  =  0,  1,  . . . ,  n-l 
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Note  that  k  can  be  written  as: 


i  +  rn  for  i  |  -  1  and  r 


•  • *  >  »  0,  1 1 


( 


so 


Exp[-j2*k^]  =  Exp[-j2i:(i+rn)^]  »  Exp[-j2iri^] 


since 


Therefore, 


Exp[j2Tt]  =  1 


2  1 


H(mw^ } 


=  E  2-.  ci+rnExp[-j2ir(i+rn)^] 
~  n  r=-» 

•  “  o 


•-1 


E  Exp[-j2iri£l  £  c. 

n  r=-»  t 

'  "  n 


+rn 


where 


£-1 


E  c'^Exp[-j2tri^] 

i  =  — 

1  O 


Ci  =  r~»  Ci+rn 


in  <  <  <  n  i 
2  -  -  2  1 


This  equation  for  H(m^)  is  the  OFT  of  the  sequence  c'..  Thus 
the  weight  coefficients  or  pulse  response  derived  by  the  frequency 


( 
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sampling  technique  are  the  c'.  which  are  "aliased"  Fourier  series 
coefficients.  QEO 


LEMMA  2 
number. 


PROOF: 


where 


f 


:  Given  a  sequence  ck  of  Fourier  coefficients  and  an 
c,  then  a  value  of  n  can  be  found  such  that  if  m  >  n 


2 


E 

|  r*  |  >0 


ci+rn 


< 


€ 


By  a  form  of  Parseval’s  theorem  [Trt] 


«  2T 

E  |ck|2  =  fj  | Hj  (tu)  1 2du  -  a 

£=-oo  -1 

2T 


H^(ui)  is  the  ideal  frequency  response 
a  is  some  constant 


arbitrary 

then 
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-  - 1 
2  2 


Z  £  ck4r„  «  E  £  |ck,r„i2 
k=f  | r } >0  k»=f  <r<>0 


<  S'si2 


k  <•—  -1 


5  I  kill  'h'Z  ' 


*  a  -  XL  |c.  |? 

\*\$ 


since  £  j c.  ( 2  is  a  convergent  series  [Ros]  there  must  exist,  for  any 
k*-«  k 

arbitrary  c,  some  n  such  that 


a  -  £  |cJZ  <  e 

w<| 


which  implies  that 


QED 
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PROOF  of  THEOREM  5-2:  the  integral  squared  error  is  defined  as: 


where 


also,  if  the  c( 

H^u)  -  Ha(t 

by  LEMMA  1 . 

Parseval ' ! 


H. (to)  is  the  ideal  transfer  function 
H  (u)  is  the  actual  transfer  function 

d 

are  the  Fourier  series  coefficients  of 


theorem  [Trt]  is: 


(w),  then 

Upt-jukT] 


where  the  ck  are  the  Fourier  series  coefficients  of  H^(u>)  so 


£  ic/  *  E 


E.  clt+rn^ 


k>|-l 


2T  *  ‘  . 

if  )H-(w)]^du>  +  53  !cifi  +  *-*  l  23  ck+rn^ 

Tj  '  1  fc,2£  lrl>0 


-1 

2T 


Given  an 


arbitrary  e  there  exists  an  n]  such  that  if  m  >  n,  then 


If  |H  U)i2du>-  23  lck>2<! 

I  '  1  -in 


km=SL 

k  2 


aiso,  by  LEMMA  2  an  n£  can  be  found  such  that  if  m  >  n?  then 


E  I  E  ‘Kt/  ‘  t 

.  -m  1 r 1 >0 
k=T 


'thus  if  n  -  MAX(nrn2)  then  for  any  m  >  n 


c2  <1*1  ’  ‘ 


I 


r 


PROOF  of  COROLLARY  to  THEOREM  5-2:  The  error  output  power  Is: 


1 

2T 


f  |E((u)  [2S(u))du> 

£1 

2T 


( 


where 


S(u>)  =  Power  Spectral  Density  of  input  signal 

S>)  =  Power  Spectral  Density  of  signal  which  is  due 
to  the  error  between  the  ideal  filter  and  the 
actual  filter 


Therefore, 

J_ 

2T 

pfi  <  /  V{£(oj)  |Zdaj  »  Ve2 

2? 

where 


V  =  Max[S(w)] 

2 

e  *  integral  squared  error  between  the  ideal  and  actual 

2 

transfer  function,  e  can  be  made  as  small  as 
desired  (Theorem  5-2) 


this  completes  the  proof. 

NOTE:  A  signal  with  zero  bandwidth,  such  as  a  sinusoid,  does  not 
have  a  finite  power  spectral  density.  Consequently  this  proof  does  not 
hold  for  sinusoids.  However,  it  is  conjectured  that  Gibbs  phenomena 
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t 


[Car]  may  appear  at  discontinuities  in  the  transfer  function.  So  i f 
the  input  has  a  sinusoid  at  the  frequency  of  a  discontinuity,  then  the 
Gibbs  phenomenon  will  cause  the  error  power  to  be  finite  for  all 
values  of  n. 

Theorem  5-2  supports  the  intuitive  notion  that  increasing  the 
number  of  weights  in  an  FIR  filter  should  allow  the  transfer  function 
to  become  closer  to  the  ideal  transfer  function.  The  corollary  lets  Its 
know  that  in  principle  frequency  sampling  filters  can  attenuate  any 
finite  bandwidth  signal  to  any  desired  level.  This  means  that  the 
filters  designed  by  the  APS3  fit  gain  tolerance  scheme  one.  (In 
practice  this  may  not  be  true  due  to  problems  such  as  finite  precision 
arithmetic). 

Figures  5-7  through  5-10  show  a  series  of  frequency  sampling 
designed  bands top  filters.  All  of  the  filters  have  a  stopband  between 
12.555  and  25%  of  the  Nyquist  frequency.  The  transition  gains,  which 
are  at  12.5%  and  25%  of  the  Nyquist  frequency,  are  0.4.  Figure  5-7 
shows  the  frequency  response  of  a  32  weight  filter.  Figures  5-8,  5-9 
and  5-10  show  the  frequency  response  of  a  64  weight  filter,  a  128  weight 
filter  and  a  256  weight  filter  respectively.  The  improvement  in  the 
frequency  response  is  evident. 

Summary  and  conclusions 

This  chapter  has  presented  a  scheme  for  an  adaptive  power  separator 
which  is  based  on  a  LALE,  a  OFT  and  an  inverse  OFT.  We  have  shown  that 
the  adaptive  power  separator,  called  the  APS3,  fits  tolerance  scheme  one. 
This  means  that  the  designer  can  select  two  power  thresholds,  and  p2 
(with  p.j  <  p2)  and  guarantee  that  all  signals  with  power  less  than 
will  be  passed  and  all  signals  with  power  greater  than  p2  will  be 


( 


stopped.  In  the  context  of  the  APS3,  and  gain  tolerance  scheme  one, 
the  signals  that  are  stopped  will  be  multiplied  by  some  gain  g  which 
is  between  zero  and  some  arbitrarily  small  value  5£  which  is  selected 
by  the  designer. 
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6  --  APS4 :  THE  ADDITION  OF  NOISE  CANCELLING  TO  THE  APS3 


This  chapter  proposes  and  examines  an  adaptive  power  separator 
(APS)  which  uses  an  adaptive  noise  canceller  to  improve  the  performance 
of  the  APS3.  The  new  APS  has  more  components  than  the  APS3,  but  the 
performance  fits  tolerance  scheme  two.  Because  of  the  significant 
improvement  in  performance,  the  new  APS  is  called  the  APS4  to  distin¬ 
guish  it  from  the  APS3. 

The  previously  explained  APS,  the  APS3,  designs  a  bandstop  filter 
that  will  reject  the  strong  signals.  However,  the  bandstop  filter  is 
not  perfect,  so  a  small  fraction  of  the  strong  si-gnal  may  still  be  in 
the  output  of  the  APS3.  Consequently,  if  the  strong  signal  is  powerful 
enough  the  output  power  can  be  very  large.  We  would  like  to  find  a 
scheme  for  cancelling  the  residues  of  strong  signals. 

Noise  cancelling  is  a  common  application  for  adaptive  filters 
[Wid,  Glo].  The  usual  situation  for  noise  cancelling  is  shown  in 
Figure  6-1.  The  adaptive  filter  forms  a  best  fit  between  the  noise 
signal  at  the  reference  input  and  the  noise  in  the  primary  input.  Thus 
the  output  consists  of  the  signal  and  whatever  noise  the  filter  could 
not  match.  In  the  case  of  the  APS4,  the  primary  input  to  the  canceller 
is  the  output  of  the  APS3,  and  the  reference  input  to  the  canceller  is 
the  output  of  a  bandpass  filter  which  passes  only  the  strong  signals. 

Figure  6-2  shows  the  structure  of  the  APS4.  The  LALE  and  the 
bandstop  filter  form  the  APS3.  The  new  bandpass  FIR  filter,  which  has 
pass  regions  wherever  the  bandstop  filter  has  rejection  regions 
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igure  6-2.  Structure  of  the  APS4. 


( 


(see  Fig.  6-3),  passes  only  the  strong  signals.  The  strong  signals  are 
used  as  the  reference  input  to  a  leaky  LMS  canceller  which  cancels  the 
residues  of  the  strong  signals  from  the  APS3  output.  (An  LUIS  canceller 
is  used  because  greatly  attenuated  versions  of  the  weak  signals  will 
pass  through  the  bandpass  filter  (due  to  ripple),  and  an  LMS  filter 
would  use  these  to  cancel  the  weak  signals  in  the  APS3  output.)  The 
output  of  the  noise  canceller  is  the  output  of  the  APS4. 

The  bandpass  filter  used  in  the  APS4  is  the  complement  of  the 
bandstop  filter  which  is  used  in  the  APS3.  A  simple  way  to  design  the 
bandpass  filter  is  to  use  the  frequency  sampling  technique.  The  ideal 
transfer  function  is  derived  from  the  ideal  transfer  fucntion  of  the 
bandstop  filter  by  replacing  all  gains  of  zero  by  gains  of  one  and  all 
gains  of  one  by  gains  of  zero.  The  resulting  bandpass  filter  (see 
Fig.  6-3)  will  pass  all  of  the  signals  which  the  bandstop  filter 
rejects  --  these  are  the  strong  signals. 

The  analysis  of  the  noise  cancelling  stage  divides  into  two  cases. 
The  first  case  is  the  analysis  of  weak  signals  and  the  second  case  is 
the  analysis  of  strong  signals. 

If  a  signal  is  weak  then  the  bandstop  filter  will  pass  the  signal 
with  a  gain  of  nearly  one,  and  the  bandpass  filter  will  pass  the  signal 
with  a  gain  of  nearly  zero.  We  will  denote  the  gain  of  the  bandstop 
filter  at  the  signal's  frequency  as  1-f^  and  the  gain  of  the  bandpass 
filter  at  the  signal's  frequency  g1 .  Figure  6-4  shows  the  canceller  in 
^  this  situation.  The  signal  power  at  the  ’’error"  output  is: 

power  Z  a2  +  (l-f1-g^g)2ps  (6-1) 
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strong  weak, 
signal  signal 


strong 
si gnal 


Nyquist 


where 


g  is  the  gain  of  the  LLMS  filter  at  the  weak 
signal's  frequency 

2 

a  is  the  noise  power  in  the  APS  input  (we  ignore 
any  power  reduction  due  to  notches  and  ripples 
in  the  bandstop  filter,  also  we  ignore  the  noise 
passing  through  the  bandpass  filter) 

Ps  is  the  power  of  the  weak  signal 

The  LLMS  algorithm  acts  as  if  there  were  white  noise  with  power  y2 
added  to  the  reference  input.  Thus  the  LLMS  filter  minimizes  the 
effective  "error"  power  which  is  given  by: 

2  2  2  2 

effective  error  power  =  a  +  (l-f-j-g^g)  p&  +  Y29  ^  (6-2) 


where 


n2 


Y29 


2  2 
n2 


the  number  of  weights  in  the  LLMS  filter 

the  noise  power  which  the  LLMS  filter  believes 
is  passing  through  it 


y2  =  a  parameter  of  the  LALE 

=  power  of  "algorithmic  noise" 

\ 

The  value  of  g  which  minimizes  the  effective  "error"  power  is: 


g* 


P-f,)g,ESNR 
g^ESNR  ♦  X 


(6-3) 
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where 


£5 

r2 


( 


I 

I 

l‘ 

I 


ESNR  = 


=  effective  signal  to  noise  ratio 

The  gain  for  the  weak  signal  from  the  primary  input  to  the  "error" 
output  is: 


0-*l)  -  9i 9* 


2  ? 

gf ESNR  +  — 
31 


(6-4) 


This- is  the  gain  which  will  apply  to  a  weak  signal  in  the  APS4. 

If  a  signal  is  strong,  then  the  bandstop  filter  will  reject  the 
signal  with  a  gain  of  nearly  zero,  and  the  bandpass  filter  will  pass 
the  signal  with  a  gain  of  nearly  one.  We  will  denote  the  gain  of  the 
bandstop  filter  as  f^  and  the  gain  of  the  bandpass  filter  as  1  -  g2- 
The  analysis  parallels  the  analysis  in  the  weak  signal  case  and  we 
find  that  the  gain  which  applies  to  the  strong  signal  is: 


(l-g2)2ESNR  +  ~ 


(6-5) 


The  reason  for  adding  the  canceller  to  the  APS3  was  to  limit  the 
output  power  due  to  strong  input  signals.  Theorem  6  shows  that  this 
is  possible. 
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2  7 

THEOREM  6:  If  y2  <  pt(l-g2)  -y,  then  the  output  power  due  to  a  signal 
with  power  p  2.  Pt  will  be  less  than  or  equal  to  the  output  power  due  to 
a  signal  with  power  p,.,  where  p  is  the  threshold  power  of  the  APS3,  g_ 
is  the  passband  ripple  in  the  bandpass  filter  and  n2  is  the  length  of 
the  LLMS  filter. 


PROOF:  The  gain  of  the  APS4  is: 


9 


I  M 

2  n. 


0-g2) 


2  JL  +  JL 


for  signals  whose  power,  p,  is  greater  than  or  equal  to  the  threshold 
power  pt- 

The  output  po wer  due  to  strong  signals  is  therefore: 


^out 


P 


The  slope  of  the  output  power  with  respect  to  the  input  power  is: 


* 


Since  p,  n2  and  y2  are  greater  than  zero,  the  slope  of  the  power-out 
versus  power-in  curve  will  be  negative  if: 


C 
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r- 

L 


=>  p  > 


2  y2 

n20-g  2r 


Thus  the  slope  of  the  pov/er-out  versus  power-in  curve  is  negative  if: 


P  > 


n20-g2r 


This  is  the  condition  that  makes  the  theorem  true.  For  if  the  power- 
out  versus  power-in  curve  has  a  negative  slope  then  an  increase  in 
input  power  will  cause  a  decrease  in  output  power.  Thus  the  theorem 
requires  that  the  slope  of  the  power-out  versus  power-in  curve  be 
negative  for  all  powers  p  _>  p^..  This  implies  that: 

Y, 


pt>* 
Kt  —  n 


20-g2r 


2  n2 


y2  <  PtO-g2>  ~2 


This  demonstrates  that  the  condition: 


2  n2 


Yj  1  Pt0'92)  2 

is  sufficient  to  guarantee  that  any  signal  with  power  p  >  Pt  will  havs 
less  output  power  than  a  signal  with  power  pt> 
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QED 


As  the  power  of  a  weak  signal  increases  the  gain  of  the  APS4 
decreases.  The  gain  variation  for  weak  signals  is  minimized  if  Yj 

is  maximized.  Since  the  APS4  should  have  a  minimum  gain  variation  for 
weak  signals,  the  best  choice  of  y2  *s: 

Yj  *  Pt(l-92)2^  (6-6> 


since  this  is  the  largest  possible  value  for  y2  that  will  still 
guarantee,  by  Theorem  6,  that  the  output  power  due  to  strong  signals 
will  decrease  as  the  input  power  increases.  By  using  this  value  for 
we  find  that  the  minimum  gain  for  weak  signals  is: 


which  implies  that 


where 


P1 


9i  +  0-g2)2 

9i  +  0-g2)z 


(6-8) 


p.j  is  the  gain  deviation  in  the  weak  signal  region 

(see  Eq.  3-2,.  the  definition  of  tolerance  scheme  two) 
g^  is  the  stopband  ripple  of  the  bandpass  filter 

g^  is  the  passband  ripple  of  the  bandpass  filter 

fj  is  the  passband  ripple  of  the  bandstop  filter 
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Since  g^ ,  and  can  be  made  as  small  as  desired  (Theorem  5-2), 

Oj  can  be  made  as  small  as  desired,  thus  the  APS4  meets  the  gain  speci¬ 
fications  of  tolerance  scheme  2  (Eq.  3-2)  in  the  weak  signal  region. 

The  maximum  gain  for  strong  signals  is: 


2  =  p2 


(6-9) 


where 


2 

P£p£  is  the  maximum  allowed  output  power  due  to 
a  strong  signal  (see  the  definition  of 
tolerance  scheme  two.  Chapter  3) 

f 2  is  the  stopband  ripple  for  the  bandstop  filter 

Since  f 2  can  be  made  as  small  as  desired  (Theorem  5-2),  p2  can  be  made 
as  small  as  desired,  thus  the  APS4  meets  the  gain  specifications  of 
tolerance  scheme  2  in  the  strong  signal  region. 

The  switching  characteristics  of  the  APS4  are  controlled  by 
the  APS3  part  of  the  APS4.  Since  the  switching  characteristic  of  the 
APS3  with  hysteresis  devices  fits  tolerance  scheme  two,  the  switching 
characteristics  of  the  APS4  also  fits  tolerance  scheme  two.  Since  the 
gain  characteristics  of  the  APS4  filters  fit  tolerance  scheme  two,  the 
APS4  fits  tolerance  scheme  two.  Thus,  any  desired  quality  of  adaptive 
power  separator,  as  measured  by  tolerance  scheme  two,  can  be  implement¬ 
ed  wfth  the  APS4. 


Summary  and  Conclusions 


This  chapter  has  presented  the  APS4,  which  consists  of  the  APS3 
with  an  adaptive  noise  cancelling  stage.  The  APS4  uses  more  computa¬ 
tions  than  the  APS3,  but  the  gain  characteristic  of  the  APS4  fits 
tolerance  scheme  two.  Because  the  APS4  controls  output  power,  rather 
than  gain,  it  is  suited  to  applications  where  the  power  of  the  strong 
signals  can  be  many  times  larger  than  the  threshold  power.  For  this 
reason,  the  APS4  is  the  best  of  the  adaptive  power  separators 
discussed  in  this  work. 
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7  --  SYNTHESIS  OF  AH  APS 

In  this  chapter  we  show  how  to  design  an  APS  given  a  problem 
specification.  The  problem  to  be  solved  is:  design  a  system  that  will 
reject  jamming  signals  and  pass  communication  signals  given  a  received 
signal  that  consists  of  narrowband  communications  signals  with  powers 
between  0.1  and  2.0,  narrowband  jammers  with  power  greater  than  3.0, 
and  white  noise  with  power  0.5.  Conventional  filters  cannot  be  used 
because  the  frequencies  of  the  various  signals  are  unknown,  and  the 
LALE  cannot  be  used  because  it  cannot  suppress  the  jammers  without  also 
suppressing  the  communication  signals. 

Because  the  problem  definition  is  incomplete,  we  have  to  make 
several  assumptions.  First,  we  assume  that  the  jammer  powers  are  less 
than  20  and  design  the  APS3  accordingly.  Second,  we  assume  that  the 
jammer  powers  at  the  output  of  the  APS  should  be  less  than  0.01.  Third, 
we  assume  that  there  will  be  at  most  20  signals  at  any  time  and  that 
the  frequencies  of  the  signals  are  between  2  MHz  and  8  MHz.  Finally, 
we  assume  that  the  jammer  powers  can  be  as  large  as  20  000  and  show  how 
to  design  an  APS4  to  accommodate  this  case. 

The  design  of  the  adaptive  power  separator  consists  of  a  number 
of  steps: 

1)  Selection  of  n,  the  number  of  LALE  weights. 

Since  an  n  weight  LALE  can  handle  up  to  n/2  sinusoids,  n  must  be 
greater  than  40.  Also,  the  FFT,  which  is  an  efficient  implementa¬ 
tion  of  the  OFT,  requires  that  n  be  an  integral  power  of  2,  so  n 
is  chosen  to  be  64. 


2)  Selection  of  T,  the  sampling  period. 


t 

The  LALE  will  give  good  spectral  estimates  only  if  the  frequencies 
of  all  the  signals  are  between  20%  and  80%  of  the  Nyquist 
frequency.  Thus  T  has  to  be  50nS  MHz).  The  sampling  rate,  in 
conjunction  with  the  number  of  LALE  weights,  determines  the  bin 
width  of  the  LALE  (313  kHz).  Any  signal  with  bandwidth  less  than 
a  tenth  the  bin  width  can  be  reasonably  approximated  as  a  sinusoid, 
so  the  APS  will  treat  all  signals  with  bandwidth  less  than  30  kHz 
as  narrowband  signals  and  behave  as  predicted.  Signals  with  wide 
bandwidth  cause  unpredictable  operation  of  the  APS. 

3)  Selection  of  y,  the  "algorithmic  noise"  power. 

The  gain  of  the  LALE  is  related  to  y,  the  sum  of  the  true  noise 

power  and  the  signal  power.  To  prevent  the  spectral  sidelobes  of  ( 

v 

powerful  signals  from  causing  errors  requires  that  a  large  LALE 
gain. correspond  to  a  signal  at  the  threshold  power.  However,  for 
good  switching  performance  the  gain  should  have  a  high  sensitivity 
to  power  changes,  which  implies  that  a  small  LALE  gain  correspond 
to  the  threshold  power.  Accordingly,  a  compromise  is  reached  by 
choosing  a  LALE  gain  of  0.5  to  correspond  to  the  threshold  power. 

For  simplicity  the  threshold  power  is  taken  to  be  pt  *  2.5  (the 
average  of  the  maximum  signal  power  and  the  minimum  jammer  power). 

Using  this  information  we  can  solve  for  y: 


C 
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0.5 


gain  at  the  threshold  power 


( 


2  2 
Y  +  a 

„  P* 

1  + 


2,2 
y  +  o 


32 


2.5 

y  +  0.5 


1  +  32 


2.5 

Y  +  0.5 


=>  Y  =  79-5  *.  (7-1) 

4)  Selection  of  n  ",  the  length  of  the  FIR  filter. 

Selection  of  n'  controls  the  maximum  gain  in  the  stopband  of  the 
FIR  filter  and  also  the  quality  of  the  power  spectral  estimate. 

The  APS3  uses  the  LALE  gain  as  a  power  spectral  estimate.  However, 
the  APS3  uses  a  finite  length  OFT  to  determine  the  gain  of  the 
LALE.  The  outputs  of  the  OFT  are  samples  of  the  LALE  gain,  and 
these  samples  may  not  coincide  with  the  peaks  in  the  LALE  gain 
that  are  caused  by  signals.  If  the  OFT  samples  bracket  a  peak, 
then  the  APS3  will  act  as  if  there  are  two  smaller  signals.  This 
effect  is  a  combination  of  leakage  and  the  “picket- fence"  effect 
[Ber].  The  picket-fence  effect  is  caused  by  the  gain  characteris¬ 
tics  of  the  bank  of  bandpass  filters  that  constitute  the  DFT 
(see  Chapter  4,  Eq.  4-2).  If  the  signal  is  on-bin,  the  bandpass 
filter  passes  it  with  a  gain  of  one,  but  if  the  signal  is  between 
bins  then  the  gain  of  the  bandpass  filter  is  0.64.  The  LALE  gain 
can  be  sampled  more  frequently  if  the  LALE  weight  vector  is  padded 
with  n'-n  zeros  and  an  n'-point  DFT  is  used  to  evaluate  the 
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gain  [Bor].  This  technique  of  "zero- fill  ing"  reduces  the  picket- 

fence  effect.  The  minimum  gain  for  an  n'  point  DFT  of  or.  n  weight 

LALE  is  S  y-X~  where  x  =  .  Thus  if  a  signal  is  present  that 

^  n 

causes  the  LALE  gain  to  peak  with  a  value  g,  as  n'  point  DFT  will 

Si  n  ( pr~r) 

return  a  peak  value  between  g  and  2n' - g.  To  find  tne 

s  nm  3 

minimum  allowable  DFT  gain  we  first  determine  that  the  peak  gain 
of  the  LALE  for  a  communication  signal  with  power  2.0  (largest 
power  to  keep)  is  0.444,  and  that  the  peak  gain  for  a  jammer  of 
power  3.0  (smallest  power  to  discard)  is  0.545.  Since 
0.444  =  0.80  x  0.545,  the  minimum  gain  of  the  DFT  must  be  greater 
than  0.80  to  avoid  confusing  weak  jammers  with  strong  communication 
signals.  From  Table  7-1  we  select  n'  =  4r.  because  the  extra  gain 
gives  us  room  for  a  large  deadband  in  the  hysteresis  threshold 
device  (which  will  reduce  the  error  rate). 

The  other  constraint  was  that  the  FIR  filter  have  a  low 
enough  maximum  stopband  gain.  It  turns  out  that  a  4n  FIR  filter 
has  a  maximum  stopband  gain  of  about  -35  d8,  and  the  requirement 
is  for  a  maximum  stopband  gain  of  =  5  x  10-^  or  -33  dB.  Thus 
the  4n  filter  fits  all  of  the  constraints,  so  n'  is  chosen  to  be 
256. 

5)  Selection  of  and  a^,  the  low  threshold  and  high  threshold  of 
«.  the  hysteresis  device. 

Once  the  minimum  DFT  gain  has  been  determined  it  is  possible  to 
find  the  minimum  OFT  output  amplitude  that  corresponds  to  a 
jammer: 


133 


Mi n [a j ]  =  Min[DFT  gain]  x  LALE  gain  for  weakest  jammer 

=  0.531 

Also,  the  maximum  DFT  output  due  to  a  communication  signal  is: 

Max[a$]  *  0.444 

In  the  absence  of  other  information  it  seems  reasonable  to  select 
the  two  threshold  values  as  1/4  of  the  difference  between 
MinC'a^]  and  Max[a's]  (see  Fig.  7-1): 

a^  =  0.466 
ah  =  0.509 

6)  Selection  of  p. 

From  Eq.  5-34  the  probability  of  error  for  the  threshold  device 


P{error}  * 


VAR{SLflLE} 


Min{ah-Max[as],  MinCa^j-a^}* 


(7-2) 


From  Eq.  5-24,  VAR{S. <_  2npa  ,  so  Eq.  7-2  becomes: 


P(error)  <  7  600p 


(7-3) 


If  we  want  P(error)  <  ]%,  then: 


V  <_  1 .3  x  10* 


(7-4) 


This  value  of  p  also  guarantees  the  stability  of  the  adaptive 
process. 
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Another  constraint  on  y  is  the  expected  design  cycle  time. 

Using  Eq.  5-33,  we  can  find  E{Tc)  >  2  500,  which  is  acceptable. 

If  the  expected  design  cycle  time  was  unacceptable,  y  would  have 

to  be  reduced. 

This  completes  the  design  of  the  APS3.  The  design  parameters  for  the 
APS3  are  summarized  in  Table  7-2. 

Now  assume  that  the  jammer  powers  can  range  from  3  to  20  000.  In 
this  case  the  APS3  will  not  be  sufficient  to  reject  the  jammers  because 
the  stopband  ripple  is  too  large.  To  obtain  the  desired  performance 
we  will  add  a  noise  cancelling  stage  to  the  APS3  and  thus  design  an 
APS4.  The  APS3  portion  of  the  APS4  will  have  the  same  parameters  as  in 
Table  7-2,  except  that  y  will  be  1 .5  x  10“^  to  guarantee  that  the  small 
v  assumptions  used  in  the  analysis  of  the  LALE  remain  valid  (this  y  is 
a  tenth  of  the  stability  limit).  The  LLMS  filter  used  for  noise  cancell¬ 
ing  will  have  as  many  weights  as  the  LALE  (i.e.  64),  and  a  y  of 
1.5  x  10" 7.  is  computed  from  Eq.  6-6  to  be  96. 

This  completes  the  design  of  the  APS4. 

Simulations  of  the  APS3  and  APS4 

Both  the  APS3  and  the  APS4  were  simulated  on  a  digital  computer. 
Figures  7-2  through  7-5  show  characteristic  curves  of  both  power 
separators  for  a  single  sinusoidal  input.  Figure  7-2  shows  output- 
power  versus  input-power  as  the  input  power  ranges  from  0.1  to  10  000, 
and  Fig.  7-3  gives  an  expanded  view  of  the  transition  region.  Figure 
7-4  shows  gain  versus  input-power  as  the  input  power  ranges  between 
0.1  and  10  000,  and  Fig.  7-5  gives  an  expanded  view  of  the  transition 
region. 
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number  of  LALE  weights 

n 

= 

64 

equivalent  LALE  noise  power 

Y 

= 

79.5 

adaption  constant 

y 

= 

1.3  x  10 

low  amplitude  threshold 

a£ 

= 

0.466 

high  amplitude  threshold 

an 

= 

0.509 

number  of  FIR  Filter  weights 

n' 

= 

256 

samp! i ng  frequency 

fs 

33 

20  MHz. 

Table  7-2.  Parameters  of  the  APS3. 
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Figure  7-3.  Expanded  output  vs.  input  power  curves  for  APS3  and  APS4. 
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Figures  7-2  and  7-3  show  how  well  the  APS4  controls  the  output 
power.  Also,  the  sharp  threshold  is  evident,  as  well  as  the  fact  that 
the  APS4  easily  meets  the  design  goals.  From  Fig.  7-3  it  is  possible 
to  determine  that  the  APS3  meets  the  design  goals  over  the  limited 
power  range  that  was  assumed  for  the  design  of  the  APS3.  (The  frequency 
of  the  sinusoid  used  to  test  the  power  separators  was  chosen  to  give 
the  poorest  possible  results;  at  other  frequencies  the  output  power  and 
filter  gain  would  have  been  lower  for  powerful  signals.) 

For  comparison.  Fig.  7-6  shows  a  gain  versus  input-power  curve  for 
a  LALE  power  separator  (APS1).  The  LALE  was  designed  to  have  a  threshold 
power  of  2.5  which  gives  the  best  possible  signal  to  jammer  enhancement. 
This  figure  shows  the  very  slow  transition  from  high  gain  to  low  gain 
that  is  characteristic  of  the  LALE. 

Figures  7-7A  and  7-7B  show  the  results  of  simulating  the  APS3  with 
many  input  sinusoids.  In  this  example  there  are  10  input  sinusoids. 

Table  7-3  lists  the  power  and  frequency  of  these  sinusoids.  Figure  7-7A 
shows  a  theoretical  power  spectrum  of  the  input  signal  and  Fig.  7-78 
shows  a  theoretical  power  spectrum  of  the  output  signal.  These  power 
spectra  are  computed  from  the  power  of  the  true  input  noise,  the  power 
of  the  sinusoids  and  the  gain  of  the  APS3.  Figure  7-7B  shows  that  the 
APS3  has  deleted  the  signals  with  power  3  or  more,  and  has  left 
unaffected  the  signals  with  power  2  or  less.  Table  7-3  also  lists  the 
output  powers  at  the  signal  frequencies  to  show  the  exact  effects  of  the 
vAPS3. 

Simulations  have  also  been  used  to  determine  the  frequency 
resolution  of  the  APS3.  It  has  been  determined  that  a  jammer  that  is 
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Frequency 
(3»  of  Nyquist) 


Input  power 


12.5 

16.5 

20.5 
28.0 
32.0 
36.0 
40.0 
44.0 
48.0 


2.0 

3.0 

2.0 

3.0 

2.0 

3.0 

2.0 

3.0 


l 
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Output  power 

2.00 
0.00065 
1 .98 
0.00029 
2.00 
0.00014 
2.06 
O.OG098 


2.0 


2.04 


within  one  "bin  width"  (i.e.  a  frequency  of  of  a  signal  can  cause 
the  APS3  to  erroneously  delete  the  signal.  This  is  caused  by  a  side- 
lobe  of  the  power  spectral  estimate  of  the  jammer  adding  to  the  spectral 
estimate  of  the  weak  signal  and  causing  the  total  estimate  to  exceed  the 
threshold  value.  Thus  for  good  performance  of  the  APS3,  the  signals 
should  be  separated  by  at  least  1.5  bin  widths  (4.9%  of  the  flyquist 
frequency  for  a  64  weight  APS). 
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8  —  SUMMARY  AND  CONCLUSIONS 


This  thesis  has  developed  the  concept  of  an  adaptive  power 
separator  and  analyzed  three  implementations  of  a  generic  power 
separator  structure.  A  good  example  of  an  APS  is  the  APS3.  The  APS3 
uses  a.LALE  to  analyze  the  input  signal,  and  a  OFT  of  the  LALE’s  weight 
vector  to  produce  a  power  spectral  estimate.  .The  outputs  of  the  DFT- 
will  be  large  at  the  frequencies  where  the  input  signals  are  strong, 
and  small  elsewhere.  The  OFT  outputs  are  passed  through  threshold 
devices  that  have  zero  output  whenever  their  input  exceeds  a  predeter¬ 
mined  threshold  and  unit  output  otherwise.  The  outputs  of  the  threshold 
devices  form  a  desired,  or  ideal,  frequency  response.  A  filter  with 
this  frequency  response  would  pass  weak  signals  with  unit  gain  and 
reject  strong  signals  with  zero  gain.  This  ideal  frequency  response  is 
converted  into  a  weight  vector  for  a  finite  impulse  response  filter  by 
a  design  technique  called  frequency  sampling.  The  frequency  character¬ 
istics  of  the  resulting  filter  can  be  made  as  nearly  ideal  as  required 
by  appropriately  selecting  the  length  of  the  filter. 

Since  the  LALE-DFT  combination  gives  good  spectral  estimates,  and 
since  the  finite  impulse  response  filter  can  be  made  nearly  ideal,  the 
APS3  can  be  designed  to  have  nearly  ideal  input-to-output  characteris- 
\  tics.  However,  there  are  other  important  measures  of  APS  performance 
such  as  response  time  or  frequency  resolution.  Also,  the  analysis 
assumed  that  the  input  to  the  APS  consisted  of  a  sum  of  narrowband 
signals  in  white  noise.  This  may  accurately  model  many  situations. 


148 


i 

It 


such  as  rejection  of  CW  jammers  or  interference  from  rotating  machinery, 
but  simulations  show  that  the  APS3  and  APS4  fail  to  function  when  the 
input  consists  of  a  wideband  signal  in  white  noise.  Thus  the  APS3  and 
APS4  would  provide  no  protection  against  wideband  jammers. 

A  possible  topic  for  future  work  would  be  the  design  of  an  APS 
that  could  reject  wideband  signals  as  well  as  narrowband  signals.  Such 
an  APS  might  use  one  of  the  modern  spectral  estimation  schemes,  such  as 
maximum  likelihood  or  maximum  entropy,  instead  of  the  LALE-OFT  combina¬ 
tion  that  is  used  in  the  APS3  and  APS4. 
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